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Abstract. In this paper we present deep VLT mufti wavelength imaging observations of the field around the radio 
galaxy PKS f 138— 262 aimed at detecting and studying a (potential) proto-cluster centered at this radio source. 
PKS f f38— 262 is a massive galaxy at z = 2.16, located in a dense environment as indicated by optical, X-ray and 
radio observations. We had already found an over-density of Lya emitting galaxies in this field, consistent with a 
proto-cluster structure associated with the radio galaxy. In addition, we find 40 candidate Ha emitters that have 
nominal rest frame equivalent width > 25 A within 1.8 Mpc and 2000 km s _1 of the radio galaxy. Furthermore, we 
find 44 objects with I — K > 4.3. This number of extremely red objects (EROs) is about twice the number found 
in blank field ERO surveys, suggesting that some EROs in this field are part of the proto-cluster. The density 
of Ha emitters and extremely red objects increases towards the radio galaxy, indicating a physical association. 
From comparisons with other K band, ERO, Ha and Lya surveys, we conclude that PKS 1138—262 is located in 
a density peak which will evolve into a cluster of galaxies. 
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1. Introduction 

The search for clusters at high redshift has two main in- 
centives: distant clusters can be used to constrain cosmo- 
logical models and they provide a reservoir of high redshift 
galaxies, which can be used to study galaxy formation and 
evolution. 

According to hierarchical clustering theories, clusters 
form by the gravitational amplification of primordial den- 
sity fluctuations. In a low density universe, fluc tuations 
cease to grow after a redshift z ~ (f^ 1 - 1) l|Peeblesl 
resulting in a cluster population that evolves very 
slowly at low redshift. In an O = 1 universe density fluc- 
tuations continue to grow even at the present epoch, im- 
plying t hat the cluster p opulation would still be evolving 
rapidly ijEke et al.lll996|h The detection of even a single 
distant massive cluster, such as MS 1054—03 at z = 0.83, 
constrains the parameters (lln, <Jh) of cosmo logical models 
jBahcall fc Fanlll998tlPonahue et alJll998|) . 

The study of galaxies in nearby and distant clusters 
provides strong constraints on their evolution and for- 
mation. It has been shown that clusters at high redshift 
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can contain a di fferent galaxy popu l ation mix than those 
at low redshift. iButcher ~k, Oemlerl (|l984|) found that at 
0.1 < z < 0.5, compact clusters have significant num- 
bers of blue galaxies, the fraction increasing with redshift, 
while at z < 0.1 cluster cores are essentially devoid of 
these. This effect can be explained if we assume a process 
in which spirals lose their gas and ability to form stars 
or are converted into early type galaxies in the course of 
their evolution. Massive ellipticals, however, dominating 
cluster cores up to z ~ 1 see m to form a homogeneous 
population. Even at z = 1.27. Ivan Dokkum et al T il200lf) 
observe a scatter of the colour-magnitude relation of clus- 
ter galaxies similar to that in lower redshift clusters. 

The tight colour-magnitude relation observed out to 
high redshifts ca n be explained by seve ral galaxy for- 
mation scenarios. iDe Propris et alJ l)l999(l propose a pas- 
sive luminosity evolution model, where galaxies form all 
their stars in a single burst at z = 2 after a monolithic 
collap se (e.g. lEggen et alJ Il96^ . Ivan Dokkum k. Franxl 
ll200ll) show that the observations are also consistent 
with a scenario in which early-type galaxies are contin- 
uously transformed from spiral galaxies, causing the pro- 
genitors of the youngest, low redshift early-type galaxies 
to drop out of the sample at higher redshifts (progeni- 
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tor bias). Models bv Ivan Dokkum fc Franxl l|200l|) show 
that about half of the early-type galaxies were morpho- 
logically transformed at z < 1 and their progenitors may 
have had roughly constant star formation rates before 
transformation. These models indicate a mean luminos- 
ity weighted formation of stars in early-type galaxies of 
z = 2.0 (for f^M = 0.3 and f^A = 0.7), consistent with the 
currently favoured hierarchical galaxy formation models 
predicting the me rging of smaller galaxies at high redshift 
l|Kauffmannll 1 99d) . A third explanation is given by metal- 
licit y differences in brigh t and faint ellipticals. A model 
bv lKauffmann fc Chariot) {l998) which includes hierarchi- 
cal formation of ellipticals out of disc galaxies which have 
formed stars at modest rates and allows for the ejection of 
metals out of discs by supernova explosions predicts the 
establishment of a mass-metallicity relation among both 
late- and early-type galaxies. In this model, large ellipti- 
cals are more metal-rich because they are formed from the 
mergers of larger discs. 

The strongest observational constraints on these mod- 
els come from the highest redshift data. Both for the study 
of cluster and galaxy evolution a sample of clusters at high 
redshift is therefore desired. 



In recent years, much effort has been invested in the search 
for distant clusters, using both optical and X-ray observa- 
tions. At z > 0.5, it becomes difficult to identify the pro- 
jected two-dimensional over-density produced by cluster 
galaxies, because large numbers of foreground and back- 
ground galaxies reduce the density contrast in the optical 
wavelength regime. However, the J — K colour of nearly 
all galaxies out to z = 2 is a simple function of red- 
shift, because their near infrared light is dominated by 
evolv ed giant stars. The h ighest redshift cluster found to 
date llStanford et al. Ill 9971 CIG J0848+4453 at z = 1.27) 
has been discovered in a near infrared field survey as 
a high density region of objects with very red J — K 
colours. Optical spectroscopy confirmed the redshifts of 
eight members and a 4.5cr ROSAT X-ray detection con- 
firms the cluster's existence. Although the detectability of 
hot cluster gas in X-rays is sever ely reduced by cosm olog- 
ical surface brightness dimming, iRosati et all l|l998|) has 
found wit hin the ROSAT Deep Cluster Survey a cluster at 
z = 1.11 ^Stanford et al.l2002l RDCS J0910+5422) with 9 
spectroscopi cally confirmed clu ster members and a cluster 
at z = 1.26 i|Rosati et al.lll999l RX J0848. 9+4452) with 6 
cluster members confirmed. Most of the confirmed galaxies 
have red colours, consistent with passively evolved ellipti- 
cals formed at high redshift (z ~ 5). The latter cluster is 
very close to CIG J0848+4453, with which it might form 
a superstructure and is possibly in the process of merging. 
The recently started XMM Large Scale Structure Survey 
ijRefregier et al.ll20 02) . covering over 64 square degrees of 
sky, should be able to detect clusters with X-ray luminosi- 
ties of 2 x 10 44 erg s _1 out to z = 2. Despite the success 
of the near infrared and X-ray techniques, it is difficult to 
push these methods to find z ^> 1 over-densities. 



A practical way to find clusters and groups of galaxies 
at high redshift is to study fields containing luminous ra- 
dio galaxies. These can be observed up to the epoch of 
galaxy formation and efficiently selected by their steep 
spectrum in the radio regime I lottgering et alJ Il994t 
lLacv et al.lll994t iDe Breuck et al.l l200d|) . The most dis- 
tant radio galaxy found to date has a redshift of 5.2 
l|van Breugel et all Il999|) . The host galaxies of power- 
ful radio sources are amon g st the most massive a t any 
redshift (larvis et alJ l20f)ll: IDe Breuck et alJ 120021) and . 
are associated with ~ 10 9 BHs l|Lacv et all l200ll 
iMcLure fc Dunlodl2002l) . 

There has long been evidence that powerful radio 
galaxies at high redshift (HzRGs, z > 2) ar e located in 
the ce nter of (forming) clust ers of galaxies. lYates et alJ 
l|l989l) and lHill fc lIUvI l|l99ll) find that the average envi- 
ronment of 70 powerful classical double radio sources at 
0.15 < z < 0.82 is that of an Abell cluster, with some in 
environments as rich as Abell class 1. At z ~ 1 and higher, 
there is also eviden ce for galax y over-densities associated 
with radio galaxies. iBestl ( 2000) presents an analysis of the 
environments of 28 3CR radio galaxies at 0.6 < z < 1.8. 
The density of if-band galaxies in these field, their angu- 
lar cross correlation amplitude and near infrared colours 
correspond to the properties of low redshift Abell richness 
class to 1 clusters. The author concludes that many, but 
not all, powerful radio g alaxies at z ~ 1 lie in cluster envi- 
ronments. Furthermore. iNakata et al] l)200l() have applied 
a photometric redshift technique based on five optical and 
near infrared images of the field of 3C324 at z = 1.2 and 
identified 35 objects as plausible cluster members. The 
evidence extends even to z = 3.8, w here observat i ons in 
the field of radio galaxy 4C41.17 bv llvison etaH l|2000l) 
tentatively reveal a number of luminous submm galaxies 
over-dense by an order of magnitude as compared to typ- 
ical fields. 

There is also evidence that the environments of HzRGs 
are dense in terms of ambient g clS, clS expected for the 
centers of (forming) clusters. Radio c ontinuum observa- 
tions of ~ 70 radio gal axies at z ~ 2 l|Carilli et alJll997l 
IPentericci et al.ll2000bh show that 20 - 30% have large 
(> 1000 rad m -2 ) radio rotation measures (RMs). These 
RMs are most probably due to magnetized, ionized gas 
local to the radio sources and are comparable to the RMs 
of lower redshift radio galaxies which lie in t he centers of 
dense , X-ray emitting cluster atmospheres l(Tavlor et alJ 
11994) . 

Although convincing evidence for high density environ- 
ments associated with radio galaxies at z > 1 has been 
demonstrated, the redshifts of possible cluster members 
have not been confirmed by spectroscopy and it is there- 
fore impossible to provide a velocity dispersion for these 
structures. In a program which is currently being carried 
out with the VLT, we are targeting the fields of luminous 
radio galaxies at z > 2 and observe these with the aim of 
detecting line emitting galaxies in the associated cluster. 
At z > 2 the Lya line is rcdshifted into the optical wave- 
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length region, where we can use narrow band filters (~ 
1%) to isolate its flux from the sky background. We have 
selected ten luminous radio galaxies at 2.2 < z < 5.2, 
for which we will carry out both imaging and multi ob- 
ject spectroscopy. The survey is progressing very well and 
has already produced the discovery of the most distant 
struc ture of galaxies known (at z = 4.1. IVenemans et alJ 
120021) . 

We consider the Lya imaging and spectroscopy as a 
first step to the characterization of the cluster properties. 
The Lya emitters in a cluster form only a fraction of the 
galaxies present and might not be representative for some 
of the cluster or galaxy properties. The program will there- 
fore be followed up by e.g. broad band imaging in several 
colours. In this paper, we present new optical and near 
infrared observations of the field of PKS 1138—262 with 
the aim of uncovering populations of Ha emitting galaxies 
and EROs in the structure. 

The radio galaxy PKS 1138-262 at a redshift of 
2.16, was selected from a compendium of more than 150 
z > 2 radio galaxies as the optimum object for begin- 
ning a high redshift cluster search. It combines most of 
the above mentioned cluster indications with a redshift 
suitable for both Lya and Ha imaging. The magnitude 
of 1138—262 is the brightest of all known radio galaxies 
close to z = 2. After correction for possible non stellar 
components, the K band magnitude is 16.8, from w hich 
a ste llar mass of 10 12 M was inferred ijPentericci et alJ 
Il997|) . The radio galaxy possesses a giant (~ 120 kpc) and 
luminous Lya nebula, with a wealth of structure: a bright 
region associated with the radio jet and filaments ex- 
tendin g over > 40 kpc l|Pentericci et al.lll997t Ikurk et all 
l2000bj) . The optical counterpart of the radio galaxy is ex- 
tremely clumpy and resolved i nto many components by 
the HST (|Pentericci et a,lll998j) . These clumps have prop- 
erties similar to LBGs. The morphology of the system is 
consistent with hierarchical models of galaxy formation in 
which the LBG building blocks will merge into a single 
massive system, such as the massive galaxies observed at 
the centers of some rich clusters . The extremely distorted 
radio morphology l|Carilli et al.l Il997i) is strong evidence 
that the jets have been deflected from their original direc- 
tion by a dense and clumpy medium. The observed rota- 
tion measures of the radio emission ( 6200 rad m -1 , the 
largest in a sample of 70 HzRGs, see ICarilli et al.l fl997: 
IPentericci et alJlioOObh and its steep gradient over the ra- 
dio galaxy components also testify that the radio source 
is embedded in a dense magnetized medium. Additional 
evidence for a dense surrounding medium comes from 
Chandra X-ray observa tions, which reveal thermal emis- 
sion from shocked gas l|Carilli et alJ 1200*2)) . The pressure 
of this hot gas is adequate to confine the radio source. 

Narrow band imaging of redshifted Ly a emission 
of a 7 'x7' region around the radio galaxy l|Kurk et alJ 
l2000al Paper I) and subsequent Lya spectroscopy 
ijPentericci et alj|2000al Paper II) revealed 14 Lya emit- 
ting galaxies and one QSO. The galaxies have redshifts 
in the range 2.16±0.02 with a velocity dispersion substan- 



tially smaller than expected for a random sample of galax- 
ies selected by the narrow band filter. In addition, the 
Chandra X-ray observations of the field o f 1138—262 have 
revea led at least five AGN at z ~ 2.16 I Pentericci et alJ 
l2002h . On the basis of the evidence from the radio galaxy 
properties, the Lya halo and the galaxy over-density, 
we concluded that the structure of galaxies surrounding 
PKS 1138—262 is (the progenitor of) a cluster. 

The new observations of 1138—262 are reported in Sect.|3 
Detection and photometry of objects in the field of 
1138—262 are presented in Sect.[3J Subsequently, the selec- 
tion from these objects of K band galaxies, EROs, candi- 
date Ha emitters and candidate Lya emitters is presented 
in Sect. The properties of the EROs and candidates are 
analyzed in Sect. [5J A discussion of the implications of 
these results for the nature of the structure can be found 
in Sect. EJ which is followed by a summary of the results 
and conclusions in Sect. Throughout this article, we 
adopt a Hubble constant of Ho — 65 km s _1 Mpc _1 and 
a A dominated cosmology: 12 m = 0.3 and 12a = 0.7. The 
over-densities of galaxies at high redshift, which have not 
yet reached virialization and/or a colour-magnitude rela- 
tion with a red sequence, but will later form clusters, will 
be called proto- clusters here. 

2. Observations and data reduction 

2.1. Optical observations 

The observations of PKS 1138—262 were carried out with 
the VLT 1 . With the aim of detecting Lya emitting galax- 
ies at z = 2.16, we have observed the field of 1138—262 
for half an hour in B band and four hours in a 2% nar- 
row band, using FORSl at Antu (UT1). Subsequent multi 
object spectroscopy of candidate emitters was also car- 
ried out with FORSl, employing three masks with inte- 
gration times of 4 to 6 hours. The optical imaging and 
spectroscopy observations are described in detail in Paper 
I and II. For an overview of both old and new observations, 
see Tabled 

We have complemented the original optical imaging 
with broad band observations in R and /, using FORS2 
at Kueyen (UT2) in 2001. The detector of FORS2 was 
a Tektronix thinned and anti-reflection coated CCD with 
2048x2048 pixels and a scale of C'2 per pixel in stan- 
dard resolution mode, yielding a field size of ~ 6.'8x6f8. 
Six exposures of 5 minutes during non-photometric condi- 
tions were taken through the i?_Special filter, which has a 
central wavelength of 6550 A and FWHM of 1650 A. The 
-fiLSpecial filter has higher transmission than the standard 
Bessel R filter and its transmission curve is almost sym- 
metrical around the central wavelength, while the Bessel 
filter has its peak at 6000 A and declines towards the 
red. Six non-photometric exposures of 7.5 minutes were 



1 Based on observations carried out at the European 
Southern Observatory, Paranal, Chile, programmes P63.0- 
0477(A&D), P65.O-0324(B), and P66.A-0597(B&D). 
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taken in service mode through the Bessel / filter, which 
has a central wavelength of 7680 A and FWHM of 1380 A. 
During visitor time, three weeks later, an additional eigh- 
teen photometric exposures of 5 minutes were observed. 
The observations were made employing a jittering pat- 
tern with offsets < 20" between exposures to minimize flat 
fielding problems and to facilitate cosmic ray removal. The 
seeing on the resultant images and the 3<r limiting magni- 
tude in a 1" aperture as measured on the central square 
arcminute of the combined images is listed in Table El 

Image reduction was carried out using the IRAF 2 re- 
duction package. The individual frames were bias sub- 
tracted, flat fielded with twilight flats and cosmic rays 
were removed. The frames were combined using the 
DIMSUM 3 package. DIMSUM builds a cumulative sky 
frame from 6 to 10 subsequent unregistered images. 
Objects in the unregistered frames were detected with 
SExtractor (|Bertm fc Arnoutsl [l996) and masked during 
the process of background determination. The obtained 
sky frames were subtracted from the images. The image 
offsets were determined by measuring the positions of a 
number (~ 20) of stars on each frame. Pixels on the CCD 
which were significantly discrepant in each sky frame were 
marked as bad pixels and the exposure time for each pixel 
was computed by DIMSUM. This exposure map is later 
used as a weight map for object detection and photometry. 
In the last step of the process, all individual broad band 
frames were combined by averaging while identified cosmic 
rays and bad pixels were omitted. The registered narrow 
band images were combined by computing the average of 
each pixel stack and rejecting pixels whose intensity levels 
were 10tx above or below the noise level expected from the 
CCD gain and readout noise specifications. 

For the flux calibration of the p hotometric data 
the standard stars GD10 8 (|Qkd Il990|) and LTT4816 
l Hamuv et all 1 19921 Il994j) were used. The / band data 
obtained during non-photometric conditions was scaled 
to the photometric / band data and the R band data 
was calibrated using an olde r R band image of 1138—262 
from IPentericci et alJ l)l997|) . Astrometric calibration was 
carried out by identifying: 18 stars in the USNO-A2.0 cat- 
alogue llMojiet^e^LL . 1998), which is tied to the Tycho cat- 
alogue l|Hoeg et al.lll997h . The absolute astrometric accu- 
racy obtained in this way is ~ 0'.'2. 

Note that the narrow and broad band images obtained 
in 1999 (Paper I) were reduced again. This time we used 
DIMSUM and obtained a homogeneously reduced set of 
images in all observed optical and infrared bands. To over- 
come differences in geometrical distortion, all images were 
mapped to match the R band image. In this way, we have 



2 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Uni- 
versities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 

3 DIMSUM is a set of scripts to reduce dithered images con- 
tributed to IRAF and written by P. Eisenhardt, M. Dickinson, 
A. Stanford, J. Ward and F. Valdes. 



obtained a first set of images with their original spatial 
resolution. A second set was made in which the /, B and 
narrow band ([On]/8000, called NB0.38 from now on) 
images were convolved with the kernel required to match 
their point spread functions to the R band image, which 
had the worst seeing conditions during observations. The 
pixel-to-pixel alignment in the final images is accurate to 
within a pixel (0'.'2) over the entire image. 

2.2. Infrared observations 

In 2000 and 2001 Antu's infrared camera ISAAC was 
employed to carry out imaging of the field of 1138—262 
in several near infrared broad bands and a narrow band 
that included the redshifted Ha emission line (He I filter, 
NB2.07 from now on). The short wavelength camera of 
ISAAC is equipped with a Rockwell Hawaii 1024 2 pixel 
Hg:Cd:Te array which has a pixel scale of 0.147", yielding 
a field of ~ 2.'5x2f5. This is significantly smaller than the 
field covered by the optical observations, so we carried out 
two pointings in the infrared: one centered at the position 
of the radio galaxy as the optical observations (a,<5j2ooo 
= ll h 40 m 48 s , -26°29T0", hereafter Fl) and one to the 
North East (a, <5 J2 ooo = H h 40 m 57 s , -26°28'48", hereafter 
J-2) covering six confirmed Lya emitters (Paper II). 

All infrared images were taken in jitter mode, where 
the telescope is offset randomly between exposures but 
never farther from the original pointing than 20". In J s , 
H and K s individual frames were exposed for 100 to 225 
seconds using sub-integrations of length 10 to 45 seconds 
to avoid over-exposure of the background. The narrow 
band frames were exposed for 300 or 375 seconds with sub- 
integrations of 60 to 75 seconds respectively. Specifications 
(date, mode, band, integration time and pointing) of all 
observations are presented in Tabled Note that only T\ 
was observed in J s and H band. Observations in K s of T\ 
were taken in ESO period 63 (P63, 1999) and period 66 
(P66, 2001). The sensitivity in P66 had increased by 45% 
compared to P63, amongst others due to an aluminiza- 
tion of the main mirror. We have scaled the measurements 
done in P63 to P66, effectively reducing the formal expo- 
sure time in P63. The total exposure time in K s for T\ in 
terms of P66 time units is 1.6 hours. 

The infrared observations were reduced in the same 
way as the optical ones. However, the atmospheric emis- 
sion in the near infrared is variable on a time scale compa- 
rable with the exposure time of individual frames, causing 
fringing residuals in the frames after background subtrac- 
tion using the median of six to ten frames. These residuals 
had to be removed in the K s and NB 2.07 frames observed 
in 2000 by subtracting a low order polynomial fit to the 
lines and columns of the masked images. An overview of 
total exposure time, limiting magnitude and resultant see- 
ing can be found in Tabled 

The infrared images were registered with the optical 
reference image using the same pixel scale. In J-l 37 ob- 
jects were used for the alignment and in T2 40, resulting 
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Date 


M 


Tel/lnstr 


Filter 


A c 


^ f w li m 


Exp 


Time 


P 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


f2-4-f999 


V 


UTf/FORSf 


B Bessel 


429 


88 


300 


f800 


f 


f2-4-f999 


V 


UTf/FORSf 


NB0.38 


38f.4 


6.5 


f800 


9000 


f 


f3-4-f999 


V 


UTf/FORSf 


NB0.38 


38f.4 


6.5 


f800 


5400 


f 


6-3-200 f 


s 


UT2/FORS2 


R Special 


655 


f65 


300 


f800 


1 


5-3-200 f 


s 


UT2/FORS2 


/ Bessel 


768 


f38 


450 


2700 


1 


27-3-200f 


V 


UT2/FORS2 


/ Bessel 


768 


f38 


300 


5400 


f 


6-f-200f 


s 


UTf/fSAAC 


Js 


f240 


f60 


45/5 


3600 


f 


f 7-4-2000 


s 


UTf/fSAAC 


H 


f650 


300 


f3/8 


2496 


f 


3-4-f999 


s 


UTf/fSAAC 


K 3 


2f60 


270 


fO/fO 


2900 


f 


4-4-f999 


s 


UTf/fSAAC 


K s 


2f60 


270 


fO/fO 


f900 


f 


6-f-200f 


s 


UTf/fSAAC 


K a 


2f60 


270 


f3/8 


2496 


1 


f6-4-2000 


s 


UTf/fSAAC 


NB2.07 


2070 


26 


60/5 


5700 


1 


f 7-4-2000 


s 


UTf/fSAAC 


NB2.07 


2070 


26 


60/5 


6300 


f 


f9-4-2000 


s 


UTf/fSAAC 


NB2.07 


2070 


26 


60/5 


5400 


f 


f0-f-200f 


s 


UTf/fSAAC 


K s 


2f60 


270 


f3/8 


4576 


2 


ff-f-200f 


s 


UTf/fSAAC 


NB2.07 


2070 


26 


75/5 


9000 


2 


7-2-200 f 


s 


UTf/fSAAC 


NB2.07 


2070 


26 


75/5 


6000 


2 


2f-2-200f 


s 


UTf/fSAAC 


NB2.07 


2070 


26 


75/5 


2250 


2 



Notes: (f ) Date: Day - Month - Year (2) Visitor (V) or Service (S) mode (3) Telescope and instrument (4) ESO filter name (5) 
Filter central wavelength in nm (6) Filter full width at half maximum in nm (7) Exposure time for single frames in seconds; for 
IR observations DfT/NDfT where DIT is integration time for sub- integration and NDfT number of sub-integrations (8) Total 
exposure time in seconds (9) Pointing on T\ or T2 (see Sect. 12.2^ . 



Table 2. Resultant images 



Band 


T 


Size 


M Mm 


Seeing 


(1) 


(2) 


(3) 


(4) 


(5) 


NB0.38 


4.0 


46.6 


26.5 


0.75" 


B 


0.5 


46.6 


27.5 


0.70" 


R 


0.5 


46.6 


26.4 


0.85" 


I 


2.0 


46.6 


26.8 


0.65" 


Js 


1.0 


7.5 


24.8 


0.45" 


H 


0.7 


7.5 


23.8 


0.70" 


K a 


1.6 


12.5 


23.0 


0.45" 


NB2.07 


4.8 


f2.5 


22.8 


0.50" 



Notes: (f ) Broad band or narrow band (see TableQfor specifi- 
cations) (2) Total exposure time (hours) (3) Field size in square 
arcminute (4) 3<j limiting Vega magnitude in 1" aperture as 
measured on central square arcminute of image (5) Seeing on 
resultant image. 



in less than one pixel difference between all images over 
the entire field. The K s and NB2.07 images of T\ and 
T2 were merged into one rectangular mosaic image. The 
overlap in T\ and T1 gives rise to a region of about one 
square arcminute in the mosaic where the noise level is 
lowest. As a final step a second set of images was made 
matching the resolution of the reference image using six 
stars to estimate the difference in point spread function. 



3. Object detection and photometry 

3.1. Catalogue sets 

The results of the observations and reduction described 
above are a set of FORS 7' x T images in NB 0.38, B, R and 
/ and a set of fSAAC 2.'5x2.'5 images in J s , H, K s (2x) 
and NB 2.07 (2x). From these, we have created catalogs of 
detected and flux-calibrated objects in order to select (i) 
objects in K s band, (ii) EROs, (iii) candidate Ha emitters 
at z = 2.16 and (iv) candidate Lya emitters at z = 2.16. 

We have used th e SExtr actor software (v2.2.1, 
iBertin k Arnoutslll996(l for object detection and photom- 
etry. Since the background noise level varies across the 
images as a result of the dithering technique employed, 
object detection was not performed directly on the final 
reduced images, but on additional images weighted to give 
a homogeneous noise level. These were created by multi- 
plication of the reduced images by their associated expo- 
sure time maps. Only for the detection of Lya emitters, a 
homogeneous noise level image convolved to the R band 
seeing was used, where the detection sensitivity for slightly 
extended objects (C/85) is highest. A disadvantage of the 
use of convolved images for object detection is that spu- 
rious sources (e.g. remaining cosmic rays) become indis- 
tinguishable from real sources. Three regions in the Lya 
detection image that were badly affected by bright stars 
were blanked (for a total of 7.75 arcminute 2 ). The source 
extraction parameters were set so that, detected objects 
must have at least 8 connected pixels with flux in excess of 
1.5 times the background noise level of the detection im- 
age, except for the Lya detection image, where a source 
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has to have 14 connected pixels. To ensure that the colours 
are computed correctly, object photometry was done on 
the convolved images, by employing SExtractor's double 
image mode using the apertures defined on the weighted 
images. A weight map created from the square root of the 
exposure time map was used to estimate the errors in the 
photometry. 

iKronl l|l98(t) and llnfantel l|l987f) have shown that 
for stars and galaxy profiles convolved with Gaussian 
seeing, > 94% of the flux is inside the appropriately 
scaled Kron aperture. We have therefore used SExtractor's 
MAG_AUTO implementation of Kron's first moment al- 
gorithm to estimate the total magnitudes of the sources. 
The resultant magnitu des were corrected f or galactic ex- 
tinction of Ab = 0.172 l|Schlegel et al.ll998(l and assuming 
an Ry = 3.1 extinction curve, which resulted in a de- 
crease of the zero-points of B, R and I by 0.2, 0.1 and 0.1 
respectively. No changes were necessary for the infrared 
zero-points. 

SExtractor classifies the likelihood of detected objects 
to be stars or galaxies using a neural network. The resul- 
tant stellaricity index has a range from 0.0 to 1.0, where 
stars should have a value near 1.0 and galaxies a value 
near 0.0. 

To derive a list of Ha emitting candidates, object de- 
tection and aperture definition needed to be done on the 
infrared narrow band image. For our first set of cata- 
logues, the apertures were therefore defined on the un- 
convolved homogeneous noise level image associated with 
the NB2.07 image. In this way, most remaining cosmic 
rays and CCD defects are too small to be included in the 
list and we do not introduce a preference for a fixed spa- 
tial frequency. Photometry with these apertures was sub- 
sequently carried out on all eight convolved images. The 
resulting catalogue contains 479 objects, of which thir- 
teen were either spurious or not suitable for the detection 
of line emitters (e.g. bright stars, a few remaining cosmic 
rays and some image boundary defects). 

We are also interested in the population of EROs in 
the field of PKS 1138-262. These objects are old ellipti- 
cal or dusty starburst galaxies at z > 1 and may also be 
present in the proto-cluster structure. EROs have extreme 
I—K colours, i.e. they are detected in K band but are very 
faint in the optical. A second set of catalogues was there- 
fore based on the unconvolved homogeneous noise level 
image associated with the K s band image. From this set, 
we derive the K s band counts and EROs. The resulting 
catalogue contains 550 objects. 

Although candidate Lya emitters were selected in 
Paper I, we have derived a new list of Lya candidates 
based on the newly reduced B and NB0.38 images using 
selection criteria consistent with the criteria for selecting 
the Ha emitters presented in this paper. For this purpose, 
a third set of catalogues was constructed, based on the 
convolved homogeneous noise level image associated with 
the NB0.38 image. This set contains 1027 sources. 



Table 3. K s band galaxy counts 





n 


N* 


Nsco 


>4 


>5 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


18.0 


19 


2.5±0.4 




1 





18.5 


17 


3.8±0.6 




6 


1 


19.0 


29 


6.2±0.7 


4.3±0.3 


9 


1 


19.5 


34 


8.9±0.8 


6.3±0.4 


19 


2 


20.0 


58 


13.5±1.0 


8.4±0.4 


26 


5 


20.5 


69 


19.0±1.2 


11.3±0.5 


37 


7 


21.0 


59 


23.8±1.4 


16.1±0.6 


47 


11 


21.5 


84 


30.5±1.6 




59 


16 


22.0 


55 


34.9±1.7 




64 


19 



Notes: (1) Limiting K s magnitude (2) Differential counts be- 
tween K l ^ m — 0.5 and Kg (3) Cumulative number arcmin -2 of 
galaxies brighter than limiting Kl lul , Poisso n error is indicated 
(4) Sa me as (3) from blank field survey bv lScodeggio fc Silval 
( 2000) (5) Number of objects brighter than Kl iul with I-K s > 
4 (6) Same as (5) for / - K s > 5 (See Sect.lOl. 

4. Number counts and cluster candidates 

4.1. K s number counts 

K band number counts were derived from the catalogues 
based on the K s image. Table|211ists the number of sources 
per half magnitude bin a nd the cum ulative number of 
sources per square degree. iBestl i|2000t BOO from now on) 
shows that selecting K band objects with SExtractor's 
stellaricity index below 0.8 efficiently selects galaxies as 
opposed to stars. This conclusion is based on the J — K 
colour of the detected objects which is in general bluer for 
stars. Because T1 is not imaged in J, we determined the 
galaxy counts from the total number counts (550 objects) 
by selecting only those objects with stellaricity index lower 
than 0.8 (470 objects). 

The cumulative counts were compared with observa- 
tio ns of a blank field of sub stantial size (43 arcminute 2 ) 
bv lScodeggio fc Silval l)2000|) . Although these authors also 
make a distinction between stars and galaxies based 
on SExtractor's stellaricity index (0.85), it is not clear 
whether the number counts in their tables are total counts 
or galaxies only. We assume here that galaxy counts are 
listed. We observe on average 1.5±0.1 times the number 
of objects expected from the blank field survey, as illus- 
trated by the observed surface density of galaxies brighter 
than K s = 20 of 13.5±1.0 ar cmin -2 , compared with 
8.4±0.4 arcmin -2 determined bv lScodeggio fc Silval Less 
deep, but us i ng a m uch larger field are the observation of 
iDaddi et al.l l)2000|) . which determine a K s galaxy count 
of 3.29±0.07 objects arcmin -2 up to K R = 18.7 (which is 
slightly higher than IScodeggio fc Silval ' s value), while we 
find 4.5±0.6 galaxies arc min -2 up to this limit. From Fig. 
2 in IDaddi et alJ l|2000j) . it is clear that there exist con- 
siderable scatter in the observations by different authors, 
which they interpreted as due to cosmological field-to-ficld 
variations of up to a factor two. Our K s band counts are 
near the upper limit of the observed variations. 
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Fig. 1. Galaxy counts as a function of distance from the 
radio galaxy in circular bins of 20" width. The error- 
bars represent Poissonian errors. The mean density (146 
xlOOO) per square degree inside the largest circle (210", 
bin not shown) is indicated by the horizontal line. 



We have analyzed the surface density of K band se- 
lected galaxies as a function of distance from the radio 
galaxy by counting the number of galaxies in circular ar- 
eas around 1138—262. Fig.^shows the number of galaxies 
per arcsecond 2 in circular bins of 20". It is clear that the 
counts show an excess within 50" or 0.45 Mpc from the 
radio galaxy. The deviation from the mean density of the 
joined first two bins is 2.2<r. 

The richness of clusters can be assessed by counting 
the number of cluster galaxies found within a radius of 
0.5 Mpc of the central galaxy with magnitudes between 
mi and mi+3, mi being the magnit ude of the centra l 
galaxy. This value, N . 5 , is defined bvlffill fc Lilly! l|l99l|) 
and based on an earlier definition bv lAbelll |l958). The K 
magnitude measured for 1138—262 is 16.1. We can correct 
this value for line emission from Ha and N II as measured 
in the NB2.07 band, by solving for line and continuum 
contributions in the broad and narrow band. Wc obtain 



c = 



AA 



bb 



AA 



rib 



(i) 



where c, b and n are the continuum flux and flux mea- 
sured in broad and narrow band and A Abb, AA n b are the 
FWHM of the broad and narrow band, respectively. The 
K magnitude in the NB 2.07 band is 15.8 and the result- 
ing magnitude of the continuum flux therefore 16.2. Note 
that this difference of 0. 1 units of magni t ude is less than 
the 0.3 units derived in IPentericci et ahl l|l997h based on 
an assumed Lya/Ha ratio of 1. In addition to line flux, 
we have to cor r ect fo r other non stellar contributions. 
IPentericci et alJ l)l997|) attribute 0.2 magnitude units to 
the point source contribution and a maximum of 0.2 units 
to the resolved contribution. Wc conclude that the con- 
tinuum K magnitude representative of the stellar popu- 
lation in the brightest cluster galaxy is close to 16.5. At 
z = 2.16, 0.5 Mpc is equivalent to 56", and the angular 
area occupied by a disc with 0.5 Mpc radius is 2.7 arcmin 2 . 



The number of galaxies between 16.5 and 19.5 in this area 
around 1138-262 is 31. From Table 1 in BOO, we find that 
blank fields contain ~ 2.06 x 10 4 galaxies per square de- 
gree between magnitude 16.5 and 19.5 (references to the 
blank field data can be found in the caption of Fig. 6 
in BOO), or 15 galaxies within a circular area of 56" ra- 
dius. The net excess co unt around 1138—2 62 is therefore 
16±6. From Table 4 in iHill fc LilTl l|l99l|) we read that 
for clusters of richne s s 0, 1 , 2, the mean N .5 values de- 
termined bv iBahcalll l|l98lh are 12±3, 15±5 29±8. The 
mean value of N .5 in 3C radio galaxy fields ifYates et alJ 
1989) is 11±2, similar to the average number of 11 galax- 
ies in the 28 z ~ 1 radio galaxy fields found by BOO. The 
measured N .5 suggests that 1138—262 is located in an en- 
vironment with a density comparable to richness class or 
1 clusters. During the evolution of this structure, however, 
more galaxies might fall in, further increasing the richness 
of the cluster. 



4.2. Extremely red objects 

Since the discovery of extremely red objects (EROs, 
lElston et allll988l Il989l) . there has been considerable in- 
terest for these high redshift galaxies because their prop- 
erties can constrain models of galaxy formation and evo- 
lution. They are now generally believed to be either 
evolved ellipticals or dusty star bur sts at z > 1 an d 
have been shown to cluster strongly l|Daddi et al.ll200o() . 
Using our multi band observations, we can search for 
EROs which could form a population of galaxies asso- 
ciated with the radio galaxy at z ~ 2.16. Clusters are 
known to possess a population of elliptical galaxies which 
form a red sequen ce in a colour-magnitude diagram (e.g. 
iBower et al1ll992|) . The evolution of the elliptical galaxy 
population in clusters has been shown by numerous au- 
thors (e.g. IStanford et al.lll998h to be simple and homo- 
geneous, indicating that the stellar population that makes 
up the red sequence is formed at high redshifts (zf > 2). 
lOladders fc Yeel {2000) show that the red sequence can be 
exploited to find clusters of galaxies up to z ~ 1.4 using 
optical imaging. Basically, the cluster red sequence is as 
red as or redder than other galaxies at a given redshift and 
all lower redshifts if pr operly chosen filters st raddling the 
4000 A break are used l|Gladders fc Yedl2000() . For a clus- 
ter elliptical at z = 2.2 the 4000 A break is redshifted to 
12800 A, in the infrared J band. We have selected EROs 
based on their / — K s colour, which also targets galaxies 
at the redshift of the proto-cluster. The samples selected 
with these bands can be compared with literature data. 



4.2.1. EROs in the field of PKS 1138-262 

Figs. [21 and show plots of I — K vs K and I — K vs 
J — K for the sources with apertures defined on the K s 
band image. These plots show that there is a considerable 
number of objects with very red colours (I—K s > 4) in our 
field. An enlargement of Fig. is shown in Fig. 0] for the 
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Fig. 2. Colour-magnitude plot of I — K vs K for the 
544 sources detected on the K s band image. The sources 
within 40" of the radio galaxy are indicated by circles. 



range in I — K > 3.75. For I — K > 4.3, there are very few 
bright objects: the median K band magnitude in the range 
4.3 < I ~ K < 5.1 is 20.9 while it is 1.3 magnitude lower 
in the range 3.5 < I — K < 4.3, with the sole exception 
of the radio galaxy (7 - K = 4.8, J - K = 2.8). This 
increase in magnitude suggests that a large fraction of the 
objects redder than this limit are distant galaxies and we 
consider therefore the 44 objects with I — K > 4.3 in 
the field of 1138—262 as EROs for the remainder of this 
paper. These objects are listed in Table IA.3I However, 
since other authors use different criteria, we have listed 
the number of red objects according to several selection 
criteria in Table 

We will compare the number density of EROs in 
the field of 1138—262 with the density observed in a 
23.5 arcminute 2 area o f the Chandra Deep Field by 
IScodeggio fc Silvalll2000h . which agrees well with the den- 
sity o f I — K 8 > 4.0 objects measured by ICowie et al.l 
l|l996h . In Table the number of objects with I - K > 4 
and > 5 found in the field of 1138—262 is shown per 
K s magnitude limit. Up to K s — 21, one expects 33±4 
(6±2) EROs with I - K > 4 (5) in a blank field of 
this size, while we observe 47 (11). Note, however that 
iDaddi et alJ l|2000l) claim that the density of EROs with 
R — K s > 5 and K < 19 i n the 43 arcmin 2 CDFS, 
derived bv IScodeggio fc Silval l|2000|h is a factor of five 
sm aller than the one derived in the 700 arcmin 2 survey 
bv lDaddi et alJ |2000). This large discrepancy is not un- 
expected given the strong clustering of EROs. Although 
we are using different selection criteria and most EROs 




Fig. 3. / — K vs J — K colour-colour plot for the 320 
sources detected in the central K s band image. Non- 
detections in / or J and therefore lower limits on the 
magnitude are indicated by arrows. The radio galaxy (in- 
dicated by a diamond) has J — K = 2.8 and I — K = 
4.7. The track marks the colour of a stellar population 
at z = 2.2 formed by a 100 Myr single burst followed 
by passive evolution for ages of 0.5 Myr to 2.4 Gyr after 
the burst. Four formation redshifts of the population are 
indicated. PKS 1138—262 is indicated by a diamond. 



detected in our comparatively small 12.5 arcmin 2 field 
are fainter than K = 19, we should be careful drawing 
conclusions from this comparison. It is unclear whether 
the clustering of EROs at fainter flux levels is as strong, 
as at bright levels, but the measurements bv lDaddi et alJ 



( 2000J) show that the clustering amplitude of EROs with 
K„ < 18.5 is t wice as high as of EROs with K s < 19.2. 
Eqs. 8 and 9 in lDaddi et alJ l|2000t) prescribe the rms fluc- 
tuation of ERO counts due to cosmic variance, given the 
clustering strength of EROs. We have made the conser- 
vative assumption that K = 21 EROs are clustered as 
strongly as K = 19 EROs, resulting in an uncertainty 
on the number of EROs with I - K > 4 (5) of 11 (4). 
Therefore, we tentatively find an over-density of about a 
factor of 1.5, which might be due to a population of EROs 
in the proto-cluster at z ~ 2.2 on top of a field population 
of EROs at lower redshift. 

An important argument for the proposition that part 
of the ERO population consists of proto-cluster members 
at z = 2.2 is the gradient in the spatial distribution of 
EROs, as shown in Fig.[5J which is similar to the distribu- 
tion of K s band counts. The density within a 40" radius of 
the radio galaxy is more than four times higher than the 
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Fig. 5. Extremely red object (I — K > 4.3) counts as 
a function of distance from the radio galaxy. The EROs 
are counted in circular bins of 20" width. The error-bars 
represent Poissonian errors. The mean density (8 xlOOO) 
per square degree outside 60" (0.5 Mpc) is indicated by 
the horizontal line. 



Fig. 4. This close-up of Fig. shows the extremely red 
objects. The horizontal line denotes our ERO selection 
criterion I — K > 4.3. Error-bars have been reduced in 
size by a factor two to increase the readability of the plot. 
The solid line, diamond and arrows are described in the 
caption of Fig. 01 Three formation redshifts for the stellar 
population are indicated. 

Table 4. Extremely Red Objects counts near 
PKS 1138-262 



I — K 


# 


#c 


(1) 


(2) 


(3) 


4.3 - 5.3 


36 


27 


> 5.3 


8 


6 


4.5 - 5.5 


25 


22 


> 5.5 


7 


5 



Notes: (1) I — A' colour (2) Number of EROs in 12.5 arcminute 2 
field covered by both K band images (3) Number of EROs in 
7.1 arcminute 2 field covered by central K a band image only. 



mean density outside a 60" (0.5 Mpc) radius. This mean 
surface density of 8 x 10 3 galaxies per square degree would 
imply a number of 28 EROs in the field, roughly consistent 
with the number expecte d from the blank field observed by 
IScodeggio fc Silval ifeoOOj) . In Fig.H the objects within 40" 
from the radio galaxy are indicated by circles. It is clear 
that many of the reddest objects lie near the radio galaxy: 
there are 17 EROs within a 40" radius of 1138-262 (12 
arcmin -2 ) and 27 outside this radius (2.4 arcmin~ 2 ). The 
over-density of EROs in the field of 1138—262 is therefore 
mostly due to the red galaxies near the radio galaxy, which 
is consistent with the observed excess of EROs being due 
to a cluster population associated with the radio galaxy. 



4.2.2. EROs in the proto-cluster at z = 2.2 

What colours do we expect for galaxies at z = 2.2? We 
have computed evolutionary tracks for several stellar pop- 
ulation using the Galaxy Iso chrone Synthesis Spectra l 
Evolution Library (GISSEL93, frWual fc CharloiJIl flflfl 
The IMF we have used in the models is a lSalpeterl l|1955|) 
law with lower mass cutoff at 0.1 M Q and upper mass cut- 
off at 125 Mq. On Fig. a track is indicated for a stellar 
population at z — 2.2 formed by passive evolution after a 
100 Myr single burst model. The track starts at an age of 
5 Myr at I — K = 0.5 and leaves the plot at an age of 2.4 
Gyr (I — K — 6.5). We have also computed the colours 
for a constant star formation model, and an exponential 
star formation model with r = 1 Gyr, but these did not 
reach I — K > 4 even after 4 Gyr, the maximum age for a 
galaxy at z = 2.2 in this cosmology. The EROs in the field 
of 1138—262 have colours consistent with galaxies which 
have undergone a starburst 0.5 Gyr ago (z — 2.6) or longer 
and evolved passively thereafter. However, dust reddening 
is not taken into account in the models presented here. 
Reddening of actively star forming galaxies by dust can 
also sh ift their I — K colours into this regime as shown 
by e.g. IPozzetti fc Mannuccl ((2000) . The I-K colour of 
their elliptical template increases to ~ 7 at z > 2, while 
the colour of their dusty starburst template stays constant 
at ~ 6. We observe no objects with I — K ~ 7, but there 
are three objects with lower limits to their I — K colour 
of ~ 6, which could be z ~ 2 ellipticals and a few objects 
with I — K ~ 6. We conclude that the excess of ~ 10 - 15 
of the EROs in the field of 1138—262 is caused mostly by 
galaxies at z — 2.2 which are reddened by dust. 
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Fig. 6. Colour-magnitude diagram for 467 sources de- 
tected in the NB 2.07 image. The dot-dashed lines are 
lines of constant excess signal E (see text). Also shown 
are lines of constant EWo (for z — 2.16). Candidate Ha 
emitters are objects with EWo > 25 or 50 A and S > 2 
or 3 (see table 0. 



4.3. Candidate Ha emitters 
4.3.1. Selection procedure 

Candidate line emitting objects were selected on the basis 
of their excess narrow versus broa d ba nd flux, following 
the cr iteria of iBunker et al ] lll995h and iMoorwood et al) 
(120001) . The selected candidates fulfill two criteria: first, 
they have sufficient equivalent width (EW) and second, 
their broad band flux is significantly lower than expected 
for a flat spectrum source. Having measured the narrow 
band flux for each source, we compute the expected broad 
band flux and its standard deviation assuming a flat spec- 
trum. The error parameter E is defined as the number 
of standard deviations the measured broad band flux de- 
viates from the expe cted broad band flu x of a flat spec- 
trum source (see also lBnnker et ahlll995l) . Note that E is 
well defined for objects not detected in the broad band. 
A NB 2.07 -K s band versus NB 2.07 magnitude plot for 
the 466 bona fide objects in the Ha selection catalog (see 
Sect.|3) is shown in Fig. El Also drawn are two horizontal 
lines indicating rest frame equivalent width (EWo) of 25 
and 50 A and two curves indicating E equal to 3 and 2. 
The curves of constant E have been computed for median 
narrow band and broad band errors; the actual E of indi- 
vidual sources depends, amongst others, on the aperture 
size and local background noise. 



for clusters at high redshift 



Table 5. Properties of the samples of Ha candidates 



EWo 


E 


niNB 


# 


#H 


Hn38 


n 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


50 


3 




17 


2 


1.4±0.3 




25 


3 




23 


3 


1.8±0.4 




25 


2 




40 


3 


3.2±0.5 




50 


1 




48 


2 


3.8±0.6 




25 


1 




60 


3 


4.8±0.6 




75 


2 


19.3 


1 





0.1±0.1 


0.61±0.35 a 


100 


3 


19.5 


2 


1 


0.2±0.1 


0.12±0.05 ft 



Notes: (1) Rest frame equivalent wid th lower limit (2) S ignal 
to noise lower limit (as defined by Bu nker et alJ Il995t) (3) 
Additional selection criterion: narrow band magnitude up- 
per limit (4) Number of candidates in 12.5 arcminute 2 area 

(5) Number of candidates within the radio galaxy Lya halo 

(6) Surface density of candidates in the field of 1138—262 
(arcm in - 2 ) (7) Surface density in o ther fields: a IBunker et alJ 
ll995Tl . b Ivan der Werf et alJ feuOOft . 



We find 17 candidate emitters with rest-frame equiv- 
alent width EWo > 50 A and E > 3, all of which have 
narrow band magnitudes < 21.3. One of these objects is 
the radio galaxy, while a second is within the extent of 
the Lya halo of the radio galaxy. If we lower the selection 
criteria to EWo > 25 A and E > 2, we find 40 candi- 
dates. All of these have mNB2.07 < 21.6. In addition to 
the two objects identified above, there is one more ob- 
ject in this list within the radio galaxy Lya halo. The Ha 
halo of the radio galaxy is interesting in itself, especially 
in comparison with the Lya halo and will be studied in 
another article (Kurk et al. in preparation). Table |S] lists 
the number of candidates in both fields for several selec- 
tion criteria, while Table lA~T1 lists the K magnitude and 
emission line properties, for the candidates with EWo > 
25 A and E > 2. The NB 2.07 narrow band filter used 
also includes the [N IIJAA6548, 6584A lines at z = 2.16, 
but in what follows, we will refer to the combined Ha + 
Nil flux and equivalent width, as Ha flux and equivalent 
width respectively, unless otherwise noted. 

4.3.2. Number density of candidate Ha emitters 

In recent years the number of succe ssful searches f or Ha 
emitters in the infrared has grown. Bu nker et al ] lll995h 
have obtained near infrared narrow band images in a 4.9 
arcminute 2 field towards the quasar PHL 957, in an at- 
tempt to detect Ha emission from a damped Lya ab- 
sorber at z = 2.313. They find 3 candidate Ha emit- 
ters at 2.29 < z < 2.35 with EW > 75 A and E > 
2, brighter th an a narrow band magnitu de limit of 19.3. 
More recently, Ivan der Werf et ahl l|2000|) presented a sur- 
vey for Ha emission at redshifts from 2.1 to 2.4. They 
have observed several fields containing known damped 
Lya systems (also including the field of PHL 957) and 
radio galaxies as well as random fields for a total area 
of 55.9 arcminute 2 . They detect two radio galaxies and a 
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Fig. 7. Position plot for the 40 (17) Ha candidates with 
EW > 25 (50) A and £ > 2 (3) indicated by open (filled) 
squares. Large squares indicate candidates for which Lya 
emission has been detected (see Sect. 15. l|) . Axes are in arc- 
seconds, the radio galaxy is at the origin. The dotted boxes 
indicate the borders of the two reduced ISAAC fields. 
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Fig. 8. Candidate Ha emitter (EW > 25 A, £ > 2) 
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candidates are counted in circular bins of 20" width. The 
error-bars represent Poissonian errors. The mean density 
(7.9 xlOOO) per square degree outside 1 arcminute (0.5 
Mpc) is indicated by the horizontal line. 



damped Lya absorber in the field of PHL 957 in Ha emis- 
sion accompanied by two close emitters which had not 
been observed before. In addition, they find another can- 
didate in the field of a radio galaxy, but at a large distance 
from it (81"), which seems to be a merging galaxy. In to- 
tal, they find seven Ha emitters with EWo > 100 A and 
£ > 3 down to an area weighted narrow band magnitude 
limit of 19.5. The narrow band fluxes of these candidates 
are > 2.0 xlO -16 erg cm -2 s -1 . 

Our VLT near infrared imaging is much deeper and 
most of our candidates have fluxes below this limit. To 
compare the number density of Ha emitters in our field 
with the above mentioned authors, we have done the selec- 
tion according to their limits, including the lower narrow 
band magnitude limits imposed by their shallower obser- 
vations. We find only one and two candidates (see Table 
|5j), resulting in number densities (wit h large Poiss on er- 
rors) slightly lower and higher than iBunker et alt s and 
van der Werf et all s surveys, respectively. 

4.3.3. Spatial distribution 

Fig. [7] shows the spatial distribution of candidate Ha 
emitters. This distribution is not homogeneous over the 
observed fields. To quantify this inhomogeneity, we have 
counted the number of candidates with EWo > 25 A and 
£ > 2 in circular bins of 20" radius around the radio 
galaxy. We have taken into account the variation in sensi- 
tivity per pixel using the weight map associated with the 
NB 2.07 image. This alters the density per bin by less than 
20% in all bins compared with the unweighted computa- 
tion. The Ha candidate selection depends also weakly on 
the K s band sensitivity, but we do not expect this to have 
a large influence on the density per bin and certainly not 
on the conclusions from this plot. The density within a 
40" distance of the radio galaxy is 5.0±0.9 times higher 



than the mean density outside a 60" (0.5 Mpc) radius. 
The high number of candidates near the radio galaxy is 
consistent with 1138—262 being located in a region that 
is over-dense in Ha emitting galaxies. 

4.3.4. Contamination by other line emitters 

We do not expect that a significant fraction of our emis- 
sion line candidates arc line emitters at lower redshift, 
since there are no strong emission line redward of Ha at 
6563 A (Paa with a strength ~ 0.2 times the strength 
Ha would imply z = 0.1). It is possible that we have 
selected galaxies at z ~ 3.13 with strong [Oin]A5007 A 
emission or galaxies at z ~ 4.55 with [On]A3727 A emis- 
sion. The luminosity distance for objects at z ~ 4.6 is 2.5 
times larger than for objects at z ~ 2.16. If any of the 
emitters is indeed an [On] emitter at this redshift, its in- 
trinsic luminosity is six times higher than that assumed 
for Ha in Table IA.1I It is therefore unlikely that these 
objects contaminate our sample significantly. More likely 
is the contribution of [O in] emitters to our sample, illus- 
trated by follow-up spectroscopy to narrow band observa- 
tions similar to ours. All six galaxies suspected to be Ha 
emitters at z ~ 2.2 from low signal-to-n oise near-infrared 
spectroscopy bv lMoorwood et all l(2000^ were found to be 
[O in] emitters at higher redshift when observed at high 
enough signal-to-noise to detect the [Oiii]A4959 A line (P. 
van der Werf, private communication). Although on the 
basis of the present data, we cannot separate Ha emitters 
at z ~ 2.2 from line emitters at higher redshift, the spatial 
distribution of the candidate Ha emitters indicates that 
at least ~ 12 are associated with the radio galaxy and are 
therefore bona fide Ha emitters surrounding the cluster. 
There is, however, one candidate which stands out by its 
brightness (K = 17.7) and its size (~ 3") from the oth- 
ers (K > 19.3). The narrow band excess emission from 
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Fig. 9. Colour-magnitude diagram for 1018 sources de- 
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Fig. 10. Position plot for the 40 (11) Lya candidates with 
EW > 15 A and E > 3 (5) indicated by open (filled) 
squares. Large squares indicate candidates for which Ha 
emission has been detected (see Sect. 15. It . Axes are in arc- 
seconds, the radio galaxy is at the origin. The dotted boxes 
indicate the borders of the two reduced ISAAC fields. 



this candidate is most likely caused by Paa emission at 
1.875 /xm redshifted to z = 0.104 and marked as such in 
Table EU 

4.4. Lya candidates 
4.4.1. Selection procedure 

Using the catalogue with apertures based on the con- 
volved and partly blanked NB0.38 image, candidate Lya 
emitting objects were selected on the basis of excess nar- 
row versus broad band fluxes, similar to the method used 
to select Ha candidate emitters. Seven spurious objects 
near the edges of the images, one saturated star and one 
blended object were removed from the catalogue. The dif- 
ference between the narrow and broad band magnitude 
against the narrow band magnitude for the 1018 definite 
objects is plotted in Fig. A horizontal line indicates 
rest frame equivalent width (EWo) of 15 A and curves 
are drawn for E equal to 5 and 3. The curves of con- 
stant E have been computed for median narrow band 
and broad band errors; the actual E of individual sources 
depends, amongst others, on the aperture size and local 
background noise. We find 11 candidate emitters with rest- 
frame equivalent width EWo > 15 A and E > 5. If we 
lower the signal-to-noise criterion to E > 3, we find 40 
candidates (see Table |A~2|) . One of these objects is the ra- 
dio galaxy, while two more are also within the extent of the 
Lya halo of the radio galaxy. We have checked the num- 



ber of candidates which would be selected out of the new 
catalogues using the criteria used in Paper I (EWa > 65 
A and .Fnbo.38 > 2 x 10" 19 erg cm" 2 s^A" 1 ). We find 
70 candidates, which is consistent with the 60 candidates 
found in Paper I, given that the image on which the cur- 
rent selection is done is ~ 24% larger than the image on 
which the catalogue of Paper I was based. 

4.4.2. Spatial distribution 

Fig.llOlshows the spatial distribution of the Lya candidate 
emitters. In the south east corner of the image is a bright 
star and a nearby galaxy which inhibits the detection of 
any faint Lya emitters. Although the distribution of the 
candidates is not homogeneous, there is not a strong indi- 
cation of a density concentration within 40" of the radio 
galaxy, but at distances > 120" the density is somewhat 
below the mean, as Fig. ^2 shows. The first bin in this 
plot contains six objects: the radio galaxy, three objects 
which are part of the filamentary Lya halo and two more 
which might be associated to the halo. Most other bins are 
consistent with the mean density of 3.4 xlOOO candidates 
per degree 2 . 

4.4.3. Number density of candidate Lya emitters 

To determine the galaxy over-density near PKS 1138—262, 
we would like to compare the number density of Lya emit- 
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Fig. 11. Candidate Lya emitter (EW > 15 A, £ > 3) 
counts as a function of distance from the radio galaxy. The 
candidates are counted in circular bins of 20" width. The 
error-bars represent Poissonian errors. The mean density 
(3.4 xlOOO) per square degree is indicated by the horizon- 
tal line. 



ters found in our field to the number density of Lya emit- 
ters in a blank field. In recent years, a number of surveys 
for Lya emitting objects at high redshift have or are being 
carried out. The utility of Lya selected galaxies, represent- 
ing the faint end of the galaxy luminosity fun ction, as trac- 
ers of large scale structure was illustrated bv lFvnbo et alJ 
l)200l[) . The redshifts and position s of eight objects at z ~ 
3.0 found by their Lya emission (|Fvnbo et all fcoOO) are 
consistent with this structure being a single string span- 
ning about 5 Mpc, the first signal ob served of a filament at 
high redshift l)Moller fc Fvnboll200l|) . Lar ger surveys hav e 
bee n carried out at higher redshifts: 3.4 llHu et alJU^ QS). 
4.5 llHu et al.lll998HRhoads et aljEuB. 4.9 llOuchi et alJ 
l2002h and 5.7 l|Rhoads fc Malhotrall20ol . One of the no- 
table results of these surveys is that the observed equiv- 
alent widths of the Lya emitters indicate that they are 
young galaxies undergoing their first burst of star forma- 
tion. There are currently no observations available of Lya 
emitters at z ~ 2.2 in a blank (or any other) field, although 
such a survey is underway at the Nordic Optical Telescope 
which has the necessary but rare high UV throughput 
(e.g. iFvnbo et all 12002(1 . We will therefore compare our 
results to the observations of a known over-densit y of 
LBGs at z = 3.09 imaged in Lya l|Steidel et al.ll2000t S00 
from now on) and a shallower but much larger survey at 
z = 2.42 llStiavelli et a,l J l200lh . assuming that the results 
obtained at higher redshift are also applicable at z = 2.2. 
T his assump t ion se ems to be justified by the conclusion 
of lYan et al.l l|2002() that the galaxy luminosity function 
does not evolve si gnificantly from z ~ 3 to z ~ 6 and by 
the observation of Ou chi et al that the Lya and 

UV-continuum luminosity functions of Lya emitters show 
little evolution between z = 3.4 and z = 4.9. 

S00 have selected a sample of 72 bona fide line excess 
emitters with a narrow band magnitude limit of NBab = 
25.0 and observed frame equivalent width (EWa) > 80 A. 



This magnitude limit corresponds to NB 0.38 = 24.3. The 
number of Lya candidates in the effective 43.6 arcminute 2 
field of 1138-262 with EW A > 80 A and NB 0.38 < 24.3 is 
11. The surface density of these objects is 0.25 arcmin -2 , 
about a fourth of the value (0.96 arcmin -2 ) measured by 
S00. Taking into account the redshift range of detectable 
Lya emitters corresponding to the FWHM of the NB 0.38 
filter gives a comoving volume density of 0.0011 Mpc -3 . 
The 72 candidates detected by S00 are located in a co- 
moving volume of 21041 Mpc 3 in our cosmology, resulting 
in a volume density of 0.0034 Mpc -3 . The overdensity 
of galaxies at z — 3.09 is a factor six, consistently de- 
termined by S00 from the redshift density of LBGs as 
compared with the general LBG redshift distribution and 
from the comoving volume density of the Lya candidates 
at z = 3.09 compar ed to a blank field survey at z — 3.43 
llCowie fc Hulll99Sft . Since the comoving volume density 
of candidate Lya emitters in our field is 3.1 times smaller 
than the density found by S00, we estimate the galaxy 
overdensity in the field of PKS 1138 —262 to be a factor 
2±1. A more direct comparison with lCowie fc HiJ l|l998h 
also gives a volume overdensity of 1.6±0.7. The quoted 
errors are derived from the Poisson noise on the number 
of emitters in the three fields, but the uncertainty due to 
the difference in blank field number density at z ~ 2.2 and 
z ~ 3.1, 3.4 might be larger. 

IStiavelli et al.l l|200lj) have carried out a search for 
Lya emitters in a field of 1200 square arcminutes us- 
ing a medium band filter and found 58 candidates at 
z = 2.4. Since this sample has on average a red colour 
(< B — I >= 1.8), they conclude that the emitters con- 
tain an older stellar component and have therefore un- 
dergone their major episode of star formation at higher 
redshift. The 58 candidates have continuum subtracted 
narrow band fluxes > 2.0 x 10 -16 erg cm -2 s _1 . Taking 
into account the difference in luminosity distance, we find 
3 candidates with fluxes > 2.6 x 10 -16 erg cm -2 s -1 in 
the field of 1138—262. This amounts to a surface den- 
sity of 0.07 candidates arcmin -2 , about a f actor two more 
than the 0.048 sources arcmin -2 found bv lStiavelli et all 
The 4% filter used in their survey implies a comoving vol- 
ume of 7.6 x 10 5 Mpc 3 and a comoving volume density of 
7.6 x 10 -5 Mpc -3 , while the 8 candidates near 1138-262 
yield a comoving volume density of 3.1 x 10 -4 Mpc -3 . The 
overdensity implied by the difference in comoving volume 
density is a factor 4±2. Although we have not corrected 
our sample of candidate Lya emitters for low redshift in- 
terlopers, it is evident that the field of 1138—262 contains 
an overdensity of emitters with respect to blank field. In 
Sect. IB~^1 we discuss the overdensity of the Lya emitters 
with confirmed redshifts and estimate the mass of the 
proto-cluster implied. 



4.5. Coincidence of Lya and Ha emitters and EROs 

We have now identified candidate cluster members on the 
basis of three different criteria and we are able to find 
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Table 6. Overdensities of the three samples 
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EROs 


4.7 2.5 


3.2 2.9 


3.3 


3.3 


Ha 


4.8 2.6 


2.6 2.5 


2.8 


3.0 


Lya 


2.6 1.0 


3.1 1.5 


3.1 


1.7 



Notes: (1) Sample, either EROs, candidate Lya or Ha emitters 
(2) Density in terms of background density (3) Significance (a) 
of (2). 

out whether there is any overlap between the three pop- 
ulations in the area covered by the infrared imaging. The 
radio galaxy fulfills all requirements: it is an extremely 
red object with both Lya and Ha emission. The extended 
emission line halo also contains several objects which are 
found to be either EROs, Ha or Lya emitters. Apart from 
the radio galaxy, there are two objects classified both as 
ERO and candidate Ha emitter. One of these is located 
in the central infrared field where a J magnitude is avail- 
able and is placed wi t hin th e starburst region defined by 
IPozzetti fc Mannuccill(200(t . The Ha emitting EROs have 
I — K magnitude ~ 4.4. There are no candidate Lya emit- 
ters with I— K colour red enough to be classified as EROs. 
This is con sistent with the idea that the EROs are dusty 
starbursts l|Pev et al.l F9991 for which we do not expect 
Lya emission due to the strong extinction but at least 
some Ha emission. There are several groups of EROs, Ha 
and Lya emitters close together, as can be seen in Fig. 
1121 None of the candidate Lya or Ha emitters have suffi- 
cient infrared or optical narrow band excess emission to be 
selected as a Ha or Lya candidate, respectively. A more 
elaborate discussion of Lya /Ha ratios of the candidate 
emitters is postponed to Sect. 15 .ll 



4.6. Comparison of the spatial distributions 

It seems (Figs. El milll) that the Ha candidates and EROs 
are more concentrated towards 1138—262 than the Lya 
emitters. To determine whether these differences are sig- 
nificant, we compare the concentration of galaxies towards 
the radio galaxy by measuring the height and width of the 
density peak on top of the background population. We de- 
fine the surface density of the background population as 
the density at 60" (0.5 Mpc) < R U 38 < 200" (1.8 Mpc), 
where R1138 is the distance from the radio galaxy. The 
heigth of the peak within certain values of R1138 is given 
in Table along with its significance a. For this compu- 
tation we have excluded the radio galaxy halo objects. It 
is clear that the density peaks of the EROs and Ha can- 
didates are more significant than of the Lya candidates 
for all Rn38 < 60". The former peaks are also more pro- 
nounced at short distances (R1138 < 40"), i.e. their width 
is smaller. We conclude that there is good evidence that 
the EROs and candidate Ha emitters are more concen- 
trated towards PKS 1138—262 than the candidate Lya 
emitters. 
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Fig. 13. Lya vs Ha flux for candidate Lya emitters (filled 
symbols) and Ha emitters (open symbols). Triangles in- 
dicate upper limits. The dotted line indicates a Lya/Ha 
ratio of 10, approximately the case B prediction. 



5. Properties of the candidates 

5.1. Lya/Ha ratios of the candidate emitters 

There are no objects selected as both Lya and Ha candi- 
dates, but for some candidates we have detected emission 
in the other line below the selection criteria used in this 
work. The small overlap might seem surprising at first 
since both lines are produced by the recombination of 
neutral hydrogen. However, the strength of the Ha line 
is about a factor of ten less t han the Lya line in case B 
recombination circumstances l)Osterbrocklll989|) . Because 
the lowest Ha line flux detected is ~ 0.55 x 10~ 17 erg 



cm 2 s 1 , candidate Lya emitters should have a Lya line 



-17 



erg cm 



to have a 



flux in excess of ~ 5.5 x 10~ 
detectable Ha counterpart. 

Of the 26 candidate Lya emitters with EWo > 15 A 
and £ > 3 in the area covered by the infrared observations, 
nine have line fluxes > 5.0 x 10~ 17 erg cm~ 2 s _1 . Six of 
the latter are part of the extended Lya halo of 1138—262. 
We will discuss the other three here in more detail. For 
candidate 561 a line flux of 6.8 x 10~ 17 erg cm~ 2 s^ 1 was 
derived from the imaging observations, but spectroscopic 
observations (described in paper II) indicate a flux of 4.0 x 
10 -17 erg cm -2 s _1 . For this Lya flux level, we do not 
expect to observe Ha emission. Candidate 778 has a Lya 
line flux of 67.6 x 10 -17 erg cm~ 2 s^ 1 (outside the range of 
fluxes displayed in Fig. I13|) and coincides with an object 
detected on the NB2.07 image with a Ha flux of 8.8 x 
10 -17 erg cm" 2 s _1 , resulting in a Lya/Ha ratio of 7.7. 
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Fig. 12. A 6.'8x6.'8 I band image of the field of PKS 1138—262. Candidate Lya emitters are indicated by diamonds, 
candidate Ha emitters by squares and EROs by circles. The two regions observed in NB 2.07 and K s band are indicated 
by the boxes. 



The object has an Ha EWo of 18 A and is therefore not 
included in the list of candidate Ha emitters. The emitter 
is further described in Sect. 15.41 Candidate 441 is part of 
a chain of emitters (confirmed in Paper II) , emitting both 
Lya and Ha, having a Lya/Ha ratio of 4.7. The Lya 
emitter does not coincide with a candidate Ha emitter, 
because there is no object detected at this exact location 
on the NB 2.07 band, but Ha candidate 145 is part of the 
same chain and only 1.3" away. 



Only one candidate Lya emitter with Lya line flux 
below 5.5 x 10 -17 erg cm -2 s _1 emits detectable Ha: can- 
didate 675 has a Lya flux of 3.1 x 10~ 17 erg cm" 2 s _1 
and an Ha flux of 0.7 x 10~ 17 erg cm~ 2 s _1 , resulting in 
a ratio of 4.4 (see Fig. IT3)l . The object was not selected 
as an Ha candidate emitter as its significance (£) is only 
1.7. The lower limits of the Lya/Ha ratio computed for 
the remaining objects (see Fig. 113(1 are in the range 2-7. 

For 31 of the 40 candidate Ha emitters with EWo > 25 
A and £ > 2 we detect no Lya emission. The up- 
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per limits for the Lya/Ha ratio for these objects are 
in the range 0.03 — 1.06. Excluding three objects in the 
halo of 1138—262, there are six Ha candidates for which 
Lya emission is detected. Their ratios are in the range 
0.10-0.73 (see Fig.[T5|). 

Note that the Ha and Lya emission for respectively 
the candidate Lya and Ha emitters has been measured 
in the apertures defined on respectively the NB0.38 and 
NB 2.07 images. This can cause differences in the Lya/Ha 
ratios for these objects within a factor two. This discrep- 
ancy is, however, not large enough to explain the obvious 
difference in ratios between the two types of candidates: 
candidate Lya emitters have Lya/Ha ratios > 2, while 
candidate Ha emitters have ratio < 1. It is easily under- 
stood that objects selected by Lya emission must have 
dust free sightlines and therefore Lya/Ha ratios close to 
the case B value. There is no such selection bias towards 
low Lya/Ha ratios for objects selected by Ha emission, 
but as star formation is generally accompanied by dust 
production, it is not surprising that the candidate Ha 
emitters have low Lya/Ha ratios. 

The spatial positions of the candidates for which both 
Lya and Ha has been detected are indicated on Fig.[7|and 
Fig. EH 

5.2. Star formation rates 
5.2.1. SFR estimators 

If clouds of neutral hydrogen in or near high redshift 
galaxies absorb the integrated stellar light shortward of 
the Lyman limit and re-emit this energy in nebular lines, 
such as Lya and Ha, they provide a direct, sensitive 
probe of the young massive stellar population. Since only 
stars with lifetimes shorter than 20 Myr contribute sig- 
nificantly to the integrated ionizing flux of the galaxy, 
emission line flux is a nearly instantaneous measure of 
the star formation rate (SFR). This emission line flux can 
be computed by stellar population synthesis models, but 
is very sensitive to the initial mass function (IMF) as- 
sumed sinc e almost exclusi vely stars wi th M > 1 
contribute l|Kennicuttl fl99 Si . Assuming a lSalpeterl l)l955(l 
IMF with mass limits 0.1 and 100 M Q and solar metal- 
licity, iKennicutt] derives, u s ing th e evolutionary synthesis 
models of Kennicutt et alJ fl994), the following relation: 



SFR (Moyr -1 ) = 7.9 x 10~ 42 L Hce (erg s" 1 ). 



(2) 



Another way to measure the SFR is observing directly 
the rest frame ultra violet light, which is dominated by 
young stars. The optimal wavelength range is 1250 to 2500 
A, longward of the Lya forest but short enough that older 
stellar populations do not contribute significantly. Using 
the same stellar population as described above, for a model 
galaxy with cont inuous star formation during ~ 100 Myr, 
IKennicuttl l|l998h finds that the luminosity in this wave- 
length region scales directly with the SFR: 



This conversion is, however, dependent on the age of the 
stellar population and mode of star formation. The SFR / 
L„ ratio is higher 4 in populations younger than 100 Myr, 
up to 57% for a 9 Myr old population and lower in 100 
Myr old populations with, for e xample, an exponentially 
decreasing star formation rate. iGlazebrook et ail (^.999) 
extensively discuss the dependence of luminosity of the Ha 
line, 1500A and 2800A continuum on age and metallicity 
of the galaxy stellar population. 



5.2.2. SFRs of candidate Ha emitters 

We have computed the SFRh q from the Ha emission of 
the candidate Ha emitters assuming they are at z — 2.16, 
correcting for a 25% contribution o f the [N II] 6548+6584 
A system ((Kennicutt fc Kendll983l . For the IMF used to 
compute Eq.|31 th e UV spectrum happens to be nearly flat 
in L„ llKennicuttlll998|) . The central wavelength of our / 
band observations corresponds to 2430 A in the rest frame 
of objects at z — 2.16 and is therefore suitable to estimate 
the SFR UV of galaxies in the proto-cluster. Because we 
do not correct for possible absorption by dust, both star 
formation estimators can be considered lower limits to the 
intrinsic star formation in the galaxies. 

Excluding the radio galaxy components, a QSO (see 
Sect. I5~^|l and a low-redshift interloper (see Table lA~T|) . 
the SFRs derived from the Ha (UV) emission are in the 
range 2 - 32 (3 - 52) M Q yr" 1 . The ratio of SFR Hq / 
SFR uv is in the range 0.3 - 2.5 with a mean of 0.8±0.1 
and a disper sion of 0.5 ( s ee Fig . I14fl . 

Recently, iBuat et~ ( 2002h have investigated the star 
formation rate determined by the Ha line and the UV 
flux in a sample of nearby star forming galaxies in clus- 
ters. They find a mean ratio SFR Hq / SFR uv of 0.8 ± 0.4 
and conclude that within the error bars the two SFR es- 
timators give consistent results. There is however a large 
scatter in the sample, with two galaxies which exhibit an 
observed ratio of ~ 0.15. The mean ratio for an accom- 
panying sample of 19 starburst galaxies is ~ 2, indicat- 
ing that more dust is present in these objects. These re- 
sults were obtained with conversion factors for luminos- 
ity to SFR significantly different from Kennicutt 's val- 
ues due to a higher low mass cutoff of the IMF and the 
use of another popula tion synthesis program (Starburst99, 
iLeitherer et al.lfl9 99). but the ratio of SFR values is (co- 
incidcntally) exactly equal and their results are therefore 
comparable with ours. Assuming different IMFs, low and 
high mass cutoffs and periods since the burst for model 
stellar populations, they obtain a theoretical range of ra- 
tios of 0.66 to 1.5 for dust free galaxies. Note that, al- 
though the methods give identical values for the ratios, 
we would o btain values 50% higher using the equations in 
Sect. 4.1 of lBuat et alJ l|2002|) due to the / band sampling 
2430 A in stead of 2000 A. 



SFR (Moyr- 1 ) = 1.4 x 10^ 8 L„ (erg s^Hz" 1 ) 



(3) 



There is an error in Sect. 2.2 of IKennicuttl dl99gl) : the 
SFR/L„ ratio will be significantly higher in younger popula- 
tions (Kennicutt, private communication). 
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Fig. 14. The ratio of SFR derived from line emission to 
SFR derived from the rest frame UV continuum emission 
versus I magnitude for candidate Ha emitters (open cir- 
cles) and candidate Lya emitters (filled boxes) . Large cir- 
cles indicate emitters coinciding with EROs. The three 
Lya emitters plotted at / = 26.8 are not detected in the 
I band. Components of 1138—262, two known QSOs and 
one low redshift interloper are excluded from the plot. 



The ratios of SFRh q / SFR uv for our sample are con - 
sistent with the values determined bv lBuat et alJ (2002). 
Ratios > 1 can be explained by the presence of dust in 
the galaxies which extinguishes the UV continuum more 
effectively than the Ha emission. The three candidate Ha 
emitters which overlap with two EROs are probably the 
most dusty objects in our sample (see Fig. I14f> . The rela- 
tively strong UV emission of the Ha candidates with ratios 
< 1 can indicate that the stellar population is older than 
30 Myr, when the Ha flux drops steeply, b ut UV flux pre- 
sists up to 1 Gyr IjGlazebrook et al.lll999IK Alternatively, 
these ratios could be caused by an unconventional initial 
mass function or a contribution from direct or scattered 
AGN light to the UV emission. It is also possible that the 
youngest stars which are responsible for the Ha emission 
are more enshrouded in dust t han the stars responsib le for 
the bulk of the UV emission IjMoorwood et alJ l2000 h 

5.2.3. SFRs of candidate Lya emitters 

Assuming that the Lya emission we observe is produced in 
case B recombination circumstances and not extincted by 
dust (as indicated by the Lya/Ha ratios and lower limits 
computed in Sect. I5.1J1 . the strength of the Ha emission 
of the Lya emitters is one tenth that of the Lya emission. 
We use Eq. [21 to derive the SFR from line emission and 
Eq. from UV emission for the candidate Lya emitters. 
Excluding the radio galaxy, the QSO discovered in Paper 
II and the QSO described in Sect. 15. 41 the range of SFRL yce 
is 0.4-4.3 Moyr" 1 , while the range of SFR ratios is 0.02- 
2.65 with a mean of 0.3±0.1 and a dispersion 0.5 (see 
Fig. I14JI . There are only two emitters with a ratio > 1. 
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Fig. 15. A histogram of SFR UV for both Ha and Lya can- 
didates. The horizontal axis shows the SFR in M yr _1 . 
The bin size is set to 5 MQyr^ 1 . The radio galaxy, two 
QSOs and a low redshift interloper are excluded from this 
plot. 



These are both not detected in the I band. In general, the 
SFRLya/ SFR uv ratios of the Lya candidates are much 
lower than those of the Ha candidates. This might be 
due to the resonant nature of the Lya recombination line. 
Although we expect the Lya candidates to contain almost 
no dust, Lya photons are absorbed and reemitted many 
times by neutral hydrogen before they escape the galaxies, 
increasing the chance to be absorbed by dust enormously. 
In addition, the wavelength of Lya is more than 1000 A 
lower than the restframe wavelength of 2430 A which the / 
band samples and suffers therefore more from extinction. 
The derived values of SFRL yQ can therefore be considered 
as lower limits to the SFRs of the candidates. 

5.2.4. SFRs derived from UV emission 

To assess whether the samples of Lya and Ha candidates 
have different SFRs, we consider the SFR UV for both kind 
of emitters (Fig. I15fl . excluding components of the radio 
galaxy, the two QSOs (see Sect. I5.4f> and one low red- 
shift interloper among the Ha sample. To compute statis- 
tical properties of these samples properly taking into ac- 
count the non-detection, we have empl oyed survival statis- 
tics, u sing the methods presented in iFeigelson fc Nelsonl 
l|l985h as implemented in the ASURV (1 .1) package 
l|lsobe fc Feigelsonlll990t lLavallev et al.lll992l) . The range 
of inferred SFR UV of the Lya and Ha sample is 0.8 - 
59 Mgyr -1 with a mean of 16±3 M0yr 1 and 2.9 - 
52 M©yr _1 with a mean of 14±2 MQyr -1 , respectively. 
The distribution of SFRs derived from the / band flux 
is therefore quite comparable for the two populations. 
In addition, a Peto-Pre ntice generalized Wilcoxon test 
I Prentice fc Marek1ll979l) shows that there is a 62% prob- 
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ability that the two samples have the same underlying 
distribution. 

Integrating the SFRs of the Lya and Ha emitters, we 
obtain values of 507 and 501 Moyr -1 , respectively. The 
comoving volumes in which we have found these emitters 
are 9696 Mpc 3 (2.110 < z < 2.163, 43.6 arcmin 2 ) and 
2076 Mpc 3 (2.114 < z < 2.194, 12.5 arcmin 2 ), resulting in 
SFR densities (SFRDs) of 0.052 and 0.24 M Q yr" 1 Mpc" 3 , 
respectively. Converted to an Einstein-De Sitter cosmol- 
ogy with Hn = 5 km s _1 Mpc -1 (to compare with 



iMadau et aTTll996|) . the SFRDs of Lya and Ha emitters 
are 0.070 and 0.32 M Q yr~ 1 Mpc~ 3 , respectively, a factor 
4 to 18 more than the fiducial SFRD, as estimated by 
the mass density of metals observed today divided by the 
present age of the Universe, whi ch is close to the lower 
limit of the SFRD at z ~ 2.75 l|Madau et alb . The ob- 
served SFRDs are about a factor ten lower than the SFRD 
of the nearby Coma and Abell 1367 clusters, as derived 
from the Ha luminosities of select ed galaxies in these clus- 
ters ([Iglesias-Paramo et al.ll2002ri . Note however, that the 
derived Ha luminosities of all our 35 Ha candidates are 
in excess of 10 41 ' 5 erg s _1 , somew hat above the value of 
L* fo und for Coma (10 4123 erg s~ 1 . llglesias-Paramo et all 



l2002h . If the Coma cluster were placed at z = 2.16, we 
would not have detected more than one of its members by 
its Ha emission. 

Another measure of the star formation activity is the 
ratio of SFR in a galaxy to its mass. This ratio does not 
depend on the density of galaxies in the (proto-) clus- 
ter. We derive the stellar mass of the candidate Ha emit- 
ters at z = 2.16 from their K band magnitudes as de- 
scribed in Sect. 15.31 The ratio of star formation rate to 
mass is therefore given by the EWo: SFR / M (yr -1 ) = 
7 x 10~ 12 EW n (A). Althoug h the Sloan r' filter used by 
llglesias-Paramo et al.l l|2002i) samples a lower restframe 
wavelength range (5440A - 6755A) than our K s filter 
(6408 A - 7263A), we assume here that the flat SED of 
galaxies above the Balmer break justifies the use of the 
same formula for the Coma cluster galaxies. There is no 
bias for EWo in the Coma cluster sample, because all 
Coma galaxies with known velocities (up to r' « 16.5) 
are selected. In our sample of Ha candidates, we have 
only selected those with EWo > 25A, so we have to ex- 
clude the Coma cluster galaxies with EWo < 25A. Twelve 
Coma galaxies remain in the sample, which has a mean 
EW of 76A, while the mean EW of the non-AGN Ha 
candidates is 188A. This difference indicates that the star 
formation rate per unit mass in the proto-cluster galaxies 
is on average 2.5 times higher than in the Coma cluster 
galaxies. 



5.3. Galaxy masses 

Stellar masses for the galaxies are best estimated from 
infrared magnitudes since at that wavelength range the 
influence of rare short-lived high mass stars is minimal. 
The central wavelength of the observed K s band magni- 
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Fig. 16. A histogram of the K s derived masses for both 
Ha and Lya candidates. The horizontal axis shows the 
mass in 10 9 Mo. The bin size is 10 x 10 9 M fi . 



tudes correspond to 6840 A in the rest frame of the galax- 
ies at z = 2.16. Since this wavelength is redward of the 
4000 A break, it is suitable to estimate galaxy masses al- 
though maybe not optimal. We have computed an L/M 
ratio for galaxies at z = 2.16 o bserved through this band 
using four models of GISSEL93 l|Bruzual fc Charlotll993(l 
representing stellar populations with a constant SFR, an 
exponential SFR and a 1 and 100 Mvr burst. The IMF 
used in these models is a ISa.lneterl lll95,^ law with lower 
mass cutoff at 0.1 M Q and upper mass cutoff at 125 M Q . 
Assuming an age of 0.5 Gyr (zf = 2.55) for the stellar pop- 
ulations of the candidate emitters, we obtain M/L ratios 
in the range of 2.5 - 7.2 M Q per L Q . We will use an M/L 
ratio of 4.8 M Q / L©, which is the mean ratio of the four 
models. At 0.25 Gyr and 0.75 Gyr the mean ratios are re- 
spectively 3.3 and 5.7 M Q / L . Note that the computed 
mass to light ratio of ~ 5 in solar units is higher than 
commonly used ratios at high redshift, because we do not 
use the bolometric luminosity of the galaxies but only the 
luminosity in the restframe wavelength range sampled by 
the K s band. 

Within the area covered by the K s band images, there 
are 19 non-AGN Lya candidates. Four of these (21%) 
are not detected in K s band. Of the 35 non-AGN Ha 
candidates, also 4 (11%) are not detected in K s band. 
Taking into account the non-detections, the mean flux 
densities expressed in K magnitude of the Lya and Ha 
samples are respectively 21.2±0.3 and 20.7±0.2. We have 
converted the K magnitudes of the candidates to galaxy 
masses assuming they are at z = 2.16, resulting in a 
range of masses of 3 - 75 and 3-113 xlO 9 M© with 
means of 17±4 and 26±4 xlO 9 M Q respectively. A his- 
togram of masses of the candidate emitters is shown in 
Fig. 1161 which shows that there are only two Lya emitters 
with M > 30 xlO 9 M Q , while there are ten Ha emitters 
in this range. A Peto-Prentice generalized Wilcoxon test 
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Fig. 17. Morphology of X-ray emitting candidate Lya 
emitter 778 accompanied by faint extended Lya emis- 
sion. The gray scale represents Lya line flux (narrow mi- 
nus broad band), where white is the background level 
and black is 3 times the rms of the noise in the im- 
age. Coordinates of the central position are a, <5j2ooo = 
ll h 41 m 2!39, -26°27'44"6. 

I Prentice fc Mareklll979l) . however, shows that there is a 
20% probability that the two samples have the same un- 
derlying distribution. We conclude that there is tentative 
evidence that the candidate Ha emitters are more massive 
than the candidate Lya emitters. 

The total mass of the non-AGN Lya and Ha emitters 
derived from the detected K band magnitudes is respec- 
tively 3.1 and 8.9 xlO 11 M Q . 

5.4. X-ray properties 

The properties of X-ray point sources in the field 
of PKS 1138—262 ar e presented in a separate paper 
l|Pentericci et alJl2002|) . Here, we only mention that three 
X-ray point sources exhibit Lya emission and one Ha 
emission. One of these docs not emit sufficient Lya emis- 
sion to be included in our catalogue (EW = 3.2 A, 
S = 3.6). A second (Lya candidate 968) is indeed spectro- 
scopically confirmed to be a QSO (see Paper II). A third 
(Lya candidate 778) has a very high equivalent width 
(EWo = 311 A). In addition, it displays a faint Lya halo 
with a size of about 10", visible in the convolved image 
(Fig. 117]) . It was not included in the catalogue of Paper I 
and has not been observed spectroscopically because it is 
not covered by all jittered NB0.38 images. The presence 
of a halo, high EW, high Lya line flux (6.7 x 10~ 16 erg 
cm~ 2 s _1 ) and X-ray emission are consistent with an AGN 
nature of this object. The near infrared spectrum of Ha 
candidate 215 emitting X-rays is presented in a forthcom- 
ing paper, which confirms the AGN nature of this object. 

5.5. Conclusions on galaxy properties 

Using imaging observations in six bands, we have inves- 
tigated several galaxy populations possibly present in a 



proto-cluster at z = 2.16. From a selection on / — K s 
colour, we find several EROs. The colours of these objects 
are not consistent with the colours computed for evolved 
ellipticals at z ~ 2. However, the number density of EROs 
is high and is increasing towards the radio galaxy, indicat- 
ing that some of the red objects must be associated with 
the radio galaxy structure. The EROs could represent the 
progenitors of cluster ellipticals with some star formation 
still going on. The presence of excess emission consistent 
with Ha radiation at z = 2.16 from two of the observed 
EROs supports this idea. A second sample is formed by 
the objects selected on NB2.07 — K s colour. We believe 
that most of these objects are Ha emitting galaxies at 
z = 2.16, associated with the structure of galaxies around 
the radio galaxy. The increase of the number density of 
these objects towards the radio galaxy again supports this 
view. The high EWs of some candidate Ha emitters have 
to be explained by an AGN contribution, but most emit- 
ters should be powered by star formation with a rate of 
10 — 100 M Q yr _1 . This is also true for the sample of can- 
didate Lya emitters selected on the basis of NB0.38 — B 
colour. 

We have compared the properties of the Lya and Ha 
emitters. The density of Ha emitters is higher close to 
PKS 1138—262, their Lya/Ha ratios are lower than for 
the Lya emitters and their K s band emission and implied 
masses are higher on average. We propose the following 
scenario to explain the observed differences. A larger frac- 
tion of the Lya emitters is still being accumulated from 
the environment as compared with the Ha emitters, while 
the more massive Ha emitters have been able to retain 
more metal rich gas and dust resulting in a lower Lya /Ha 
ratio. 

6. PKS 1138-262: a rich distant cluster? 

In this section we investigate the cluster properties of the 
structure associated with radio galaxy PKS 1138—262. 

6.1. The overdensity of Lya emitters in redshift space 

In Paper II, we have presented the redshift distribution 
of the 15 confirmed Lya emitters. This distribution (see 
Fig. 3 in Paper II) does not follow the sensitivity curve of 
the narrow band filter used to select the candidate Lya 
emitters, but is narrower and centered on the redshift of 
the radio galaxy, indicating that the Lya emitters have a 
physical connection to the central radio galaxy. To quan- 
tify the significance of this deviation from randomness we 
have carried out Monte Carlo simulations of the distribu- 
tion of redshifts observed through the narrow band filter. 
We have simulated 10,000 realizations of the redshift dis- 
tribution of 14 Lya emitters. The redshift of each sim- 
ulated emitter was chosen according to the probability 
function derived from the sensitivity of the NB 0.38 filter 
and the VLT overall efficiency. The mean and standard de- 
viation of each realization is plotted in Fig. El The mean 
(z = 2.155) and standard deviation (Sz = 0.0109) of the 
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redshift distribution of the confirmed Lya emitters is indi- 
cated by a cross on this figure, which diverge respectively 
2. Oct and 2.4ct from the simulated values 5 . The measured 
distribution deviates therefore 3.1a from a random distri- 
bution, meaning that the probability that the redshifts of 
the confirmed Lya emitters are drawn from a such a dis- 
tribution is less than 0.2%. This shows that the ensemble 
of Lya emitters must have a physical connection with the 
central radio galaxy. 
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Fig. 18. The mean and standard deviations of the 10,000 
simulated redshift distributions. The cross indicates the 
position of the measured redshift distribution of the con- 
firmed Lya emitters and the contours indicate the 2 and 
3ct deviation from the combined mean and standard devi- 
ation from the simulated distributions. 



6.2. Mass and overdensity of the proto-cluster 

In Paper II the mass of the proto-cluster and volume den- 
sity of the Lya emitters was shortly discussed. A more 
thorough discussion follows here. 

We can consider two limiting cases of the dynamical 
state of the structure of galaxies: it is detached from the 
Hubble flow and completely virialized or it is no dynamical 
entity at all and the galaxies move according to the ex- 
pansion of the Universe. In the first case, the velocities of 
the galaxies are representative of the potential well of the 



5 Note that the narrow band filter was designed for the cone 
angle of FORS1 and its transmission was measured within the 
instrument. The measured deviation of the mean redshift can 
therefore not be due to the dependence of interference filter 
specifications on the incident beam. 



cluster, while in the latter case the redshifts are directly 
related to distances. 

In the first case, we can compute a virial mass as an 
estimate of the cluster mass. The group of galaxies at z ~ 
2.16 consists of 15 members, including PKS 1138—262 but 
excluding the QSO at z = 2.183, which is a 3ct outlier. The 
redshift distribution of the members is bimodal around 
the radio galaxy (see Paper II). The velocities Vi of the 
galaxies relative to the mean velocity of the group are 
computed in the following way: 

Vi = c(zi-z) / (l + z). (4) 

The velocity dispersion of the two groups of 7 members 
each (the radio galaxy is excluded in this case) is 385 
km s _1 and 2 04 km s -1 , as com puted by the Gapper 
sigma method l|Beers et al.l ll990h which is considerably 
smaller (23%) than computed by the standard method 
(square root of variance). The velocity dispersion of the 
entire sample is 901 km s _1 , as computed by the Gapper 
sigma method, which is 10% smaller than by the standard 
method. 

The virial radius R v and mass M v of N cluster mem- 
bers at positions r*j with velocity dispersion a v is computed 
as follows: 

M v = 1 a 2 v Ry. (6) 

ISmall et alJ l|l998|) show that the virial mass estimator 
works reasonable well even for a bound system with sub- 
stantial substructure. If the system is only marginally 
bound, the estimate is too high but not more than a factor 
of two. The largest uncertainties in our estimate are the 
unknown dynamical state (i.e. one group or two groups, 
bound or unbound) and the small number of available red- 
shifts. 

The virial radii of the two groups are 0.8 and 1.1 Mpc, 
about two thirds of the maximum observable radius of 1.5 
Mpc which was fixed by the size of the CCD. The virial 
masses of the groups are therefore 1.7 and 0.6 xlO 14 M . 
The mass of the two groups together (2.3 xlO 14 M Q ) is 
much smaller than the virial mass computed for the group 
of 14 galaxies as a whole (14 xlO 14 M Q ). 

In the second case, we assume that the emitters just 
detached from the Hubble flow and are now collapsing 
to a common center of gravity. In this case the redshifts 
of the emitters can be used as distance indicators with a 
small correction for their motion towards each other. We 
can compute the comoving space density of Lya emitters, 
its associated galaxy over-density and from this, derive 
a mass over-density. To compute the space density, we 
only consider the 14 confirmed emitters of Paper II in the 
redshift range 2.139 < z < 2.170. The (uncorrected) co- 
moving volume corresponding to this range is 5611 Mpc 3 . 
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The comoving volume density of the confirmed emitters 
is therefore 0.0025 Mpc~ 3 , a factor 1.4 smaller than the 
density found in for the LBG peak at z = 3.09 (S00). 
Because not all candidates were observed and some were 
too faint to produce a detectable line in their spectra, the 
resulting over-density of 4.4±1.2 with respect to the field 
population of Lya emitters is a lower limit. 

The t otal mass of t he str ucture can be estimated by 
following Isteidel et all l)l99S|) 

M = p V (1 + S m ) = 6.6 x 10 15 M Q (1 + S m ) (7) 

with 7j the mean density of the universe and S m the mass 
over-density within the volume containing the confirmed 
emitters. The mass density is related to the galaxy over- 
density <5 ga i through 

1 + bS m = C (1 + S sai ), C= l + /-/(l + ^ m ) (1/3) (8) 

where b is the bias parameter, C takes into account the 
redshift space distortio ns caused by p eculiar velocities and 
/ is the velocity factor l|Peeblesll993j) . This factor depends 
on the redshift and its value is 0.96 for z = 2.16. From the 
statistics of p eaks in the redshift distribution of LBGs, 
ISteidel et ail l|l99Sj) argue that b > 4. Taking b in the 
range 3-5, <5 m is estimated to be 0.6-1.2. This implies a 
mass of our structure of (1.1 — 1.5) x 10 16 Mq. 

The age of the universe at z — 2.155 is 3.7 Gyr, while 
the crossing time for each group is about 3 Gyr. It is there- 
fore impossible for the system around 1138—262 to be in 
a relaxed state since it takes a few crossings of the cluster 
members to virialize. Hubble flow expansion of the galax- 
ies could lead to an artificially enlarged velocity dispersion. 
The virial mass of 1.4 x 10 15 M Q must therefore be con- 
sidered an upper limit to the mass of the bound (part of 
the) system. Similarly, part of the mass computed with 
the second method is not bound to the system. Following 
the same procedure for the used volume without an over- 
density results in a mass of 6.5 x 10 15 M©. This mass is 
mostly intergalactic gas which will disappear out of the 
cluster with the Hubble flow. A more conservative esti- 
mate of the bound mass at z = 2.16 is therefore 4— 8x 10 15 
Mq , although the true bound mass of the system is prob- 
ably much lower, which does not exclude the possibility 
that more mass wil be bound at a later time. 

6.3. Conclusions on cluster properties 

Using the lists of clusters members selected in the first 
part of this paper, we have investigated the cluster prop- 
erties of the galaxies in the field of PKS 1138—262. About 
5% of the galaxies in the local universe are gathered in 
groups or clusters whose space density is larger than one 
galaxy per cubic megaparsec, about two orders of mag - 
nitude greater than the average density toresslerlll98l . 
Dense and p opulous clusters (as defined and cataloged by 
lAbelllll958l) contain on the order of 100 galaxies within 
two orders of magnitudes of the third brightest member. 
The cluster members are gravitationally bound and to 



a large extent in dynamical equilibrium (i.e. virialized). 
Because the universe is about 4 Gyr old at z = 2.2, we 
do not expect large virialized systems by that time and 
we will therefore call the over-density of objects around 
PKS 1138—262 a proto-cluster which will develop into a 
cluster. The evidence for an over-density of galaxies associ- 
ated with the radio galaxy comes from the high density of 
K s band galaxies and ERO counts as compared with field 
counts and the indication that most of this over-density is 
located in a 40" region around the radio galaxy. The den- 
sity of K band sources near the radio galaxy corresponds 
to the density of current day's clusters richer than class 
or 1. Furthermore, the density of candidate Ha emit- 
ters is similar or higher than the density found around 
other known high redshift (active) galaxies and is clearly 
highest in the 40" region around PKS 1138-262. Finally, 
the comoving volume density of candidate Lya emitters 
in the field of the radio galaxy is higher than in a blank 
field, if compared to th e number of Lya emitters in the 
field at z = 2.4 and 3.4 l|Stiavelli et alJl200ll ICowie fc HiJ 
199 8]) and in the over- density of LBGs at z = 3.09 found 
bv ISteidel et alJ {2000). From the imaging alone, we con- 
clude that PKS 1138—262 is located in an over-density 
of galaxies, while the spectroscopical observations provide 
additional confirmation. The 14 Lya emitters with red- 
shifts confirmed by spectroscopy have a redshift distribu- 
tion which shows that they are clearly associated with 
the radio galaxy. The probability that we have observed 
a rare realization of randomly distributed emitters is less 
than 0.4%. The comoving volume density of the confirmed 
emitters is already higher than in the field and is only a 
subset of the emitters in the volume. We therefore con- 
clude that PKS 1138—262 is located in a density peak 
which will most probably evolve into a cluster of galaxies. 

An indication that the proto-cluster is not virialized 
comes from the lack of extende d X-ray emission perpen- 
dicular to the radio sources axis (|Carilli et all2002|) . These 
observations imply an upper limit to the 2-10 keV lu- 
minosity of an extended relaxed cluster atmosphere of 
1.5 x 10 44 erg s _1 , less than 40% of the X-ray luminos- 
ity of the Cygnus A cluster. Given the unrelaxed state of 
the system, it is not possible to make an estimate of the 
bound mass at z = 2.16, but an upper limit of ~ 10 15 M© 
is computed for a virialized system with the properties of 
the PKS 1138-262 proto-cluster. 

Finally, we consider the galaxy population in the 
proto-cluster. As discussed in Sect. 14.21 there is a pop- 
ulation of EROs present in the cluster, most plausibly 
identified with dusty star forming galaxies. Contrary to 
what is observed in low redshift clusters, we do not ob- 
serve a sequence due to passive ellipticals in the J — K vs 
K colour-magnitude diagram, presented in Fig. ^5] Both 
results are consistent with the idea that cluster galaxies 
have evolved and were more luminous and bluer at high 
redshift (e.g. Ivan Dokkum et alJll998h . 
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Fig. 19. Colour-magnitude plot of J — K vs K for the 306 
bona fide sources detected on the K s and J band images. 
The sources within 40" of the radio galaxy are indicated 
by circles and the expected locus of the sequence of passive 
cluster ellipticals at z = 2.2 is indicated by the grey bar. 
The J — K locus is computed from the model defined in 
Sect. EH2~2l (see also Fig.©. 

7. Summary 

We have presented new near infrared and optical imaging 
of the field of radio galaxy PKS 1138-262 at z = 2.16. 
From these data and published optical imaging and spec- 
troscopy, we find several galaxy populations associated 
with the radio galaxy. 

• In the 12.5 square arcminute field covered by the K s 
band images, we detect 550 objects in K s . The number 
density of K band sources is a about 1.5 times higher 
than blank field counts and increases towards the radio 
galaxy. 

• In the same area, we find 44 objects with I—K > 4.3. 
The surface density of EROs is about a factor two higher 
than the mean density found in blank fields. Because the 
density of EROs increases towards the radio galaxy, we 
conclude that some of the EROs are galaxies at z ~ 2.2 
and associated with the radio galaxy. Most of these are 
star forming galaxies reddened by dust, although up to 
five may be passive ellipticals. 

• Also in the same area, we find 40 candidate Ha emit- 
ters. The density of the brightest emitters is comparable 
to the density in other AGN fields and, again, increases 
towards the radio galaxy. The upper limit to the Lya /Ha 
ratio of these sources is 1.1, but for most candidates the 
ratio is much smaller, meaning that these objects contain 
some dust or neutral hydrogen. Two Ha candidates are 



also identified as EROs. These must be dusty starbursts 
at z ~ 2.2. 

• In the effective 43.6 square arcminute area of the 
optical images, we find 40 candidate Lya emitters. The 
surface density of these candidates does not show an in- 
crease towards the radio galaxy. The Lya/Ha ratio of the 
26 Lya candidates in the area covered by infrared imag- 
ing is in the range 2 — 7. The sight lines through which 
we observe the Lya radiation are therefore relatively dust 
free. 

• The SFRs derived from rest frame UV flux of the 
sample of Ha and Lya candidates are comparable, but 
the galaxy masses derived from the flux around ~ 7000A 
are on average larger for the Ha candidates than for the 
Lya candidates. From the spatial distribution of the two 
samples, we deduce that the Ha emitters are more settled 
into the potential well of the proto-cluster than the Lya 
emitters. 

• From the spatial distribution, number density and 
redshift distribution of objects detected in the field of 
PKS 1138—262, we conclude that the radio galaxy is lo- 
cated in a density peak which will evolve into a cluster of 
galaxies at the present day. The bound mass at z = 2.16 
is 10 15 M at most. 

The structure around PKS 1138—262 is one of the few 
proto-clusters known at z > 2 and deserves further study. 
X-ray observations of this field have been carried out by 
ICarilli et a.lJ l)2002h and an analysis of the X-ray point 
so urces and AGN conten t of the proto-cluster is presented 
in lPentericci et all l|2002jh Infrared spectroscopy of a small 
sample of the candidate Ha emitters has been carried out 
and will be reported in Kurk et al. (in preparation). 
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Table A.l. 

and £ > 2 



Candidate Ha emitters with EWq > 25 A 
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Notes: (1) Catalog number (2) Right ascension and declination 
in J2000 coordinates (3) K band magnitude (4) Rest frame 
equivalent width (5) Significance of excess flux (6) Continuum 
subtracted narrow band flux in 10 -17 erg cm ~ 2 s _1 . Notes to 
individual objects: a Magnitude and morphology indicate that 
this object is a low redshift interloper, b Inside Lya halo of 
PKS 1138-262, c QSO (see Sect.Ol. 
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Table A. 2. Candidate Lya emitters with EWo > 15 A 
and £ > 3 



Table A. 3. Extremely red objects (I - K > 4.3) 
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Notes: 


(1) Catalog number (2) Right ascension and declination 
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in J2000 coordinates (3) B band magnitude (4) Rest frame 
equivalent width (5) Significance of excess flux (6) Continuum 



subtracted narrow band flux in 10 erg cm 



individual objects: a Inside Lya halo of PKS 1138-262, 
described in Sect. 15.41 c QSO discovered in Paper II. 



Notes to 
QSO 



Notes: (1) Catalog number (2) Right ascension and declination 
in J2000 coordinates (3) K band magnitude and error (4) I—K 
colour and error. Notes to individual objects: a PKS 1138—262. 



